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GENERAL  INTRODUCTION 


Because  of  the  high  acceleration  of  a  projectile  in  a  gun  tube,  and 
because  of  the  high  precision  that  is  sought >  gun  dynamics  requires 
consideration  of  secondary  effects  that  are  negligible  in  ordinary 
structural  analysis.  Slight  deflections  of  the  tube  are  aggravated  by 
centrifugal  action  of  the  projectile.  Vertical  oscillations  of  the  tube 
are  coupled  with  sidewise  oscillations  because  of  gyroscopic  action  of  the 
projectile.  These  effects  ax'e  studied  in  this  report. 

Section  1  is  a  development  of  the  kinematics  of  a  projectile  in  a 
concentric  moving  flexible  tube.  Aside  from  the  assumptions  that  there  is 
no  balloting,  that  the  central  axis  of  the  bore  is  inextensional ,  and  that 
plane  cross  sections  of  the  tube  remain  plane,  unstrained,  and  normal  to  the 
deflected  axis  of  the  bore,  the  basic  theory  in  Sections  1  and  2  is  exact. 
Insofar  as  kinematics  is  concerned,  dynamic  unbalance  of  the  projectile, 
the  Bourdon  effect,  axial  inertia  of  the  tubo,  and  axial  friction  of  the 
projectile  are  irrelevant. 

Although  deflections  and  twist  of  the  tube  are  small,  the  engineering 
appi'oximafcions  in  beam  theory  are  precluded .  The  Ki.rchhofv-Cieb.sch  theory 
of  large  deflections  of  thin  rods  is  a  natural  starting  point  for  the 
analysis.  This  theory  is  presented  in  A.  E.  H.  Love's  "Mathematical  Theory 
of  Elasticity,"  but  in  a  manner  that  Love  acknowledges  "is  not  without 
difficulty."  Gibbs'  vector  analysis  helps  to  clarify  arguments  of  this  kind. 
It  is  used  in  this  i*eport.  However,  with  a  view  to  computer  programming,  the 
results  ai'e  expanded  in  scalar  notation. 

The  Kirchhoff-Clcbsch  theory  treats  statics  of  a  deformed  rod.  For 
problems  of  gun  dynamics,  it  must  be  extended  to  admit  time-dependent 
deflections  and  twist  of  the  tube.  A  part  of  Section  1  deals  with  this 
problem.  On  the  basis  of  kinematics  of  the  tube,  kinematic  relations  for 
the  projectile  are  derived.  Section  2  treats  the  forces  and  moments  that, 
act  on  a  px'ojectile  in  tei'ms  of  kinematic  variables  of  the  tube. 

Section  3  ti'eats  the  motion  of  a  tapered  cantilever  tube  that  is 
actuated  by  the  projectile  and  by  prescribed  motion  at  the  bi'eech.  bonding 
of  the  tube  caused  by  the  Bourdon  effect,  axial  inertia  of  the  tube,  and 
axial  friction  of  the  projectile  are  neglected,  but  the  theoi'y  can  be 
generalized  to  include  these  effects.  The  deflections  and  twist  of  the  tube 


are  represented  as  series  of  flexural  and  torsional  modes  of  n  uniform 
cantilever  beam.  The  coefficients  in  those  series  arc  generalized 
coordinates  ot  the  tube.  They  arc  functions  of  time  that  are  determined  bv 
Lagrange's  equations.  The  theory  exhibits  gyroscopic  action  of  the 
projectile  that  causes  coupling  between  vertical  und  horizontal  oscillations 
of  the  tube.  However,  quantitative  studies  of  this  phenomenon  must  await 
numerical  computer  analysis.  A  much  simpler  problem  that  displays  the  same 
characteristics  is  treated  in  Section  4,  namely,  a  rigid  curvod  tube  that  is 
hinged  at  the  breech  so  that  it  can  swing  sideways. 


KINEMATICS  01-  A  PROJECTILE  TN  A  CONCENTRIC  FLEXIBLE  TUBE 


1.1  INTRODUCTION 

Expressions  for  the  velocity,  the  acceleration,  the  angular  velocity, 
the  kinetic  energy,  and  the  virtual  work  of  o  rigid,  spinning  projectile  in 
a  concentric  flexible  tube  are  derived  in  this  section.  Approximations  arc 
deferred  to  the  last  article  (Art.  1.14) .  The  general  theory  is  exact, 
aside  from  tho  assumptions  that  there  is  no  balloting,  that  the  axis  of  the 
bore  is  inextensional ,  and  that  piano  cross  sections  of  the  tube  remain 
piano,  unstrained,  and  normal  to  the  central  axis  of  the  bore  when  the  tube 
is  bent  and  twisted.  The  axis  of  symmetry  of  the  pi'ojectile  is  assumed  to 
be  tangent  to  the  deflected  axis  of  the  bore  at  the  location  of  the 
centroid  of  the  projectile. 

Gibbs’  vector  analysis  is  used.  A  brief  development  of  vector  analysis 
that  suffices  for  the  present  applications  is  presented  in  Appendix  C  of 
Ref.  1.  A  scalar  formulation  of  the  theory  is  presented  in  Arts.  1.11, 

1.12,  and  1.13. 

1.2  NOTATIONS 

A  bar  ovor  a  letter  denotes  a  vector. 

A  caret  over  u  letter  denotes  a  unit  vector. 

A  dot  over  a  letter  denotes  tho  derivative  with  respect  to  time  t. 
d/dt  denotes  the  total  (or  substantial)  derivative  with  respect  to  time 

(Rq.  (1.20)). 

An  asterisk  *  denotes  tho  deformed  state  of  the  tube. 

Subscripts  s  and  t  denote  partial  derivatives  with  respect  to  arc 
length  s  and  timo  t,  respectively. 

For  the  next  several  notations,  refer  to  Figures  1,  2,  and  3. 

C  is  the  unde formed  axis  of  tho  bore. 

C"  is  the  deflected  axis  of  the  bore. 

1". Dynamics  of  Rigid  Guns  with  Straight  Tubes,"  BbM-AMC  Final  Report  PAAK-U- 
80-C-0039-Task  2,  Army  Fee  ear oh  and  Development  Command,  11RL,  Aberdeen 
Proving  Ground ,  Maryland. 
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j*(  R1  nro  orthogonal  unit  vectors.  In  the  directions  of  the 
principal  normal,  the  binormal,  and  the  tangent  to  curve  C,  respectively . 
j;  is  jut.  length  on  curves  C  and  C*. 

i*,  j  * ,  It*  are  unit,  vectors  that  coincide  with  lines  in  the  deformed 
tube  which  initially  have  directions  i'.  i',  k'. 

t  is  the  unit  tangent:  vector  of  curve  C*;  Cl  “  £*) ,  Figure  1. 
fi  is  the  principal  unit  normal  of  curve  C* ,  Figuro  3. 
b  is  the  unit  binormal  of  curve  C*,  Figure  3. 

1  /It  is  the  curvature  of  curve  C*. 

1/E  is  the  tortuosity  (usually  culled  "torsion"  in  differential 
geometry)  of  curve  C*. 

Tq  is  the  tortuosity  of  Curve  C. 

w  is  a  vector,  such  that  the  vector  triad  (t*,  J*,  R*)  Issuing  from  a 
point  P*  on  C*  is  brought  parallel  to  the  orthogonal  triad  (l*  +  <$i*, 
j  *  +  6j*.  R*  dR*)  at  a  neighboring  point  P^*  on  C*  by  the  infinitesimal 
rotation  w  ds,  whore  ds  is  the  distance  P*P*  Figure  2. 

—  wmm  ^  A 

k,  ic 1  ,  t  are  components  of  w,  defined  by  w  ■  i*k  +  j*tc'  +  k*t. 
a  is  the  angle  between  vectors  j*  and  fi  at  a  point  P*  on  curve  C* , 
Figure  3, 

r(s,t)  is  the  radius  vector  from  a  designated  fixed  origin  to  a  point 
P*  on  curve  C*;  (Dr/Ds  =  t) ,  Figure  1. 

w(s,t)  is  the  angular  velocity  of  a  cross  section  of  the  tuba. 
u>.,  uu ,  are  components  of  u>,  defined  by 

1  M 


co  =  1*0)^  +  j  u>2  ¥  ^  ^3  >  (03 


£(t.)  is  the  value  of  s  locating  the  centroid  of  the  projectile  at  time 
t. 

V  -  dr/dt  is  the  velocity  of  the  centroid  of  the  projectile  (liq.  (l.liO). 
£(t)  is  the  speed  of  the  projectile  relative  to  the  tube, 
a  is  the  acceleration  of  the  centroid  of  the  projectile, 
u)  is  the  angular  velocity  (spin)  of  the  projectile  relative  to  the  tube. 
(It  is  not  the  magnitude  of  vector  w.) 


The  subscript  "a",  denoting  "axial"  io  used  rather  than  the  subscript  "t ", 
denoting  "tangential"  to  avoid  confusion  with  the  subscript  "t"  denoting 
the  partial  derivative  with  respect  to  time. 


ft  is  the  absolute  angular  velocity  of  the  projectile  (Eq.  (1.30)). 

^n’  ^b’  are  comPonents  °f  defined  by  Eq.  (1.31). 
m  is  the  mass  of  the  projectile. 

i1  is  the  moment  of  inertia  of  the  projectile  about  a  transverse  axis 
through  its  center  of  mass, 

is  the  moment  of  inertia  of  the  projectile  about  its  longitudinal 

axis . 


Tp  is  the  kinetic  energy  of  the  projectile. 

F  is  the  axial  frictional  force  on  the  projectile, 
is  the  pressure  on  the  base  of  the  projectile, 
is  the  resisting  pressure  ahead  of  the  projectile. 
A  is  the  cross-sectional  area  of  the  bore. 

M  is  the  rifling  torque. 


X  = 


rt 

u)dt 


6W  is  the  virtual  work  of  the  forces  associated  with  the  projectile. 
C**  is  a  varied  curve,  lying  infinitesimally  close  to  curve  C*. 
oo^,  10^ ,  ioa  are  components  of  u),  defined  by 


to  =  no 


n 


+  Goo, 


+  tw 


5  =  (o7 

<1  .J 


Vn>  V^,  Vg  are  the  components  of  the  velocity  V  on  the  principal 
normal,  the  binormal,  and  the  tangenc  of  curve  C*. 

an*  ab’  aa  are  co,nPonents  the  acceleration  a  on  the  principal 
normal,  the  binormal,  and  the  tangent  of  curve  C*. 

ty(s,t)  is  the  angular  displacement  of  a  cross  section  of  the  tube  in  its 
plane  (see  Eq.  (1.44)). 

x,  y,  z  are  rectangular  coordinates  attached  to  a  Galilean  reference 
frame . 

i,  j,  k  are  unit  vectors  along  the  axes  x,  y,  z. 

F  denotes  the  net  force  on  the  projectile. 

_ _  AAA 

F^,  F^,  Fa  denote  the  components  of  F  in  the  n,  b,  t  directions. 

H  denotes  the  angular  momentum  of  the  projectile  with  respect  to  its 
center  of  mass. 

H  ,  H,  ,  H  denote  the  components  of  H  in  the  n,  b,  t  directions, 
n  b  a 
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M  denotes  the  moment  about  the  center  of  mass  of  the  projectile  of  all 
the  forces  that  act  on  the  projectile. 

M^,  denote  the  components  of  M  in  the  n,  b,  t  directions. 

J  denotes  the  mass  moment  of  inertia  about  the  hinge  line  of  a  rigid 
tube  and  attached  breech. 

M  is  the  gyroscopic  couple  that  the  projectile  exerts  on  a  rigid 
curved  immovable  tube  (Figure  7) . 


X  =  co  +  'f-  . 

a  l 


u(t),  v(t)  are  the  x  and  y  components  of  displacement  of  the  tube  at 
the  breech  (Eq.  (3.13)). 

0 (t)  ,  <J>(t) ,  C(t)  are  the  x,  y,  and  z  components  of  rotation  of  the  tube 
at  the  breech.  Also,  0  is  the  angle  of  the  tangent  to  the  center  line  of  a 
rigid  curved  tube  (Figs.  6,  7,  and  8). 

p  is  the  mass  density  of  the  tube. 

I(s)  is  the  moment  of  inertia  of  a  cross  section  of  the  tube  about  a 
diameter . 

S(s)  is  the  cross-sectional  area  of  the  tube,  excluding  the  bore. 

g  is  the  acceleration  of  gravity. 

X,  is  the  length  of  the  tube. 

E  is  Young's  modulus. 

G  is  the  shear  modulus. 

a  ,  3n  are  constants  defined  by  Eqs.  (3.7)  and  (3.8),  and  by  Table  1. 

f  (s)  is  the  n'th  natural  bending  mode  of  a  uniform  elastic  cantilever 
beam  (Eq.  (3.6)) . 

ij;  (s)  is  the  n'th  natural  torsional  mode  of  a  uniform  straight  tube 
n 

that  is  fixed  at  one  end  and  free  at  the  other  (Eq.  (3.11)). 

X  (t),  Yn(t),  Zn(t)  are  coefficients  in  the  modal  expansions  of  the 
deflection  and  twist  of  the  tube  (Eq.  (3.13)).  They  are  generalized 
coordinates  of  the  tube. 

L  =  T  -  U  is  the  Lagrangian  function. 


1.3  THE  KIRCHHOFF-CLEBSCH  THEORY  OF  THIN  FLEXIBLE  RODS 

The  Kirchhoff-Clebsch  theory  of  bending  and  twisting  of  thin  rods  is 
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presented  in  References  2,  3,  and  4.  In  the  present  work,  the  rod  is  taken 
to  be  the  tube  of  a  gun.  The  undeformed  axis  C  of  the  bore  is  an  arbitrary 
curve.  The  principal  normal,  the  binormal,  and  the  tangent  to  curve  C  are 
orthogonal  unit  vectors,  denoted  respectively  by  (l ' ,  j * ,  k') .  If  the 
undeformed  tube  is  straight,  (i * ,  3',  £')  may  be  any  constant  orthogonal  unit 
vectors  such  that  k'  coincides  with  line  C.  When  the  tube  is  deformed. 


curve  C  passes  into  another  curve  C*.  Nearby  points,  P  and  ,  on  C  pass 


into  points  P*  and  P*^  on  C*.  The  lengths  PP^  and  P*P*^  are  both  taken  to 
be  ds;  i.e.,  extensionality  of  the  axis  of  the  bore  is  neglected.  Love  and 
Basset  (Refs.  2  and  3)  also  assumed  inextensionality . 

Lines  in  the  tube,  issuing  from  point  P  on  c  in  the  directions  i',  j', 

a  AAA  A 

k'  pass  into  lines  in  the  directions  i*,  j*,  k*  (Figure  11.  Vector  k*  is 
the  unit  tanger  of  C*.  Since  plane  cross  sections  of  the  tube  are  assumed 
to  remain  plane,  unstr  ned,  and  normal  to  the  centroidal  axis  C*,  the 


vectors  i*,  j*,  k*  are  mutually  perpendicular.  Love  called  straight  lines 


coinciding  with  vectors  i*,  j*,  k*  the  "principal  torsion- flexure  axes"  of 

ft  ^ 

the  rod.  Although  £*  =  t,  where 

vectors  iv  and  j*  generally  do 


‘s  the  unit  tangent  of  curve  C*,  the 
:oincide  with  the  principal  normal  n  and 


the  binormal  b  of  curve  C* . 

AAA 

The  vectors  t,  n,  b  are 
thumb,  the  forefinger,  arn. 
simultaneously  pointed 
the  right-hand  con\ 


to  be  a  right-handed  system;  i.e.,  the 
e  finger  of  the  right  hand  can  be 
actions  t,  n,  and  b.  Consequently,  with 
.e  vector  product, 


b  =  t  x  n  ,  n  ■  b  x  t  ,  t  =  nxb 


(1.1) 


With  the  approximation  ds  =  ds*,  Frenet’s  formulas  in  the  differential 


2 

A.  E.  H.  Love,  The  Mathematical  Theory  of  Elasticity ,  4th  (  Cambridge 
University  Fress,  1934,  Chap.  XVIII,  pp.  381-398. 

7 

A.  B.  Bas  et,  "On  the  Deformation  of  Thi '■*  Elastic  Wires,"  American  Journal 
of  Mathematics ,  Vol.  17 ,  1895,  pp.  28 *  ’7. 

"The  Kirchhoff-Clebsch  Theory  of  Thi.  as  tic  Rods,"  Interim  Report  BLM- 

AMC-81-2,  Contract  No.  DAAK-ll-80-C-uj39,  Army  Research  and  Development 

Command,  BRL,  Aberdeen  Proving  Ground ,  Maryland. 

u  .  a  ,  'i 

Love  considered  a  rod  of  arbitrary  cross  section,  and  he  took  i  ’  and  j  ’  to 

be  along  the  principal  axes  of  inertia  of  the  cross  section. 


geometry  of  curves  (Ref.  5)  are 


9$  n  3n 

as  "  r  *  as 


1  9£ 

r  *  as 


(1.2) 


where  1/R  is  the  curvature  and  1/E  is  t.he  tortuosity  of  curve  C*.  Partial 
derivatives  are  indicated  in  Eq.  (1.2),  because  t,  n,  and  b  generally 
depend  on  time  t  as  well  as  on  the  arc  length  s.  However,  in  this  article, 
t  is  a  passive  parameter,  since  a  single  configuration  of  the  tube  is 
considered. 

If  the  point  P*  on  C*  moves  to  a  neighboring  point  P*^  on  C*  (Figure  2), 
while  the  curve  C*  is  unchanged,  the  vectors  i*,  j*,  k*  receive  increments 
6i*,  6]*,  6ft*,  such  that  i*  +  <$!*,  j*  *  <5 j * ,  and  k*  +  6k*  are  mutually 
orthogonal  unit  vectors.  This  transformation  could  be  accomplished  by  a 
translation  and  a  rigid-body  rotation  of  the  sytem  (i*,  j*,  k*)  .  It  is 
shown  in  the  kinematics  of  a  rigid  body  that  an  infinitesimal  angular 
displacement  is  a  vector  quantity  (Ref.  6) .  Consequently,  there  is  an 
infinitesimal  vector  w  ds  which  represents  the  rotation  that  brings  the 
system  (i*,  j*,  k*)  into  parallelism  with  the  system  (i*  +  6i*,  j*  +  6)*, 
k*  +  6k*),  where  ds  is  the  distance  P*P*1  (Figure  2).  Love  defined  scalars 
K,  tc '  ,  x  by 


~W  =  i*K  +  j  *K  *  +  K*T 


(1.3) 


The  orthogonal  projection  of  curve  C*  onto  the  j*k*  plane  (or  the  k*i*  plane) 
is  a  plane  curve  v.ith  curvature  K  (or  <’)  at  point  P*.  The  tortuosity  of 
the  undeformed  axis  C  is  denoted  by  xQ.  It  is  the  rate  of  rotation  of  the 

osculating  plane  of  curve  C  with  respect  to  arc  length  s.  The  deformational 

twist  of  the  tube  per  unit  length  is  accordingly  t  -  tq. 

If  a  rigid  body  undergoes  an  infinitesimal  angular  displacement  w  ds 
about  a  fixed  axis,  and  if  p  is  a  radius  vector  from  a  point  on  that  axis  to 

a  particle  Q  of  the  body,  the  displacement  of  Q  is  w  x  p  ds.  Letting  p 

stand  successively  for  i*,  j*,  and  k*,  and  noting  Eq.  (1.3),  we  get 


5D.  Struik,  Differential  Geometry }  Addison-Wesley  Press,  Cambridge,  Mass., 
1950 , 

^E.  T.  Whittaker,  Analytical  Dynamies,  4th  ed. ,  Dover  Publications ,  New 
York,  1944,  Chap. ~Tf~Avt.  8. 
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=  w  X  1*  =  3*t  -  £*<» 

tts 


f)i  •“  A  A  A 

7^—  =  W  X  3*  =  k*K  -  1*T 


=  W  X  k*  =  i**1  -  j*K 


(1.4) 


A  A  A  A 


Figure  3  shows  the  vectors  i*,  j*,  n,  b  in  the  normal  plane  of  curve 
C*.  It  follows  from  the  geometry  of  the  figure  that 


i*  =  n  sin  a  -  b  cos  a  ,  3  *  =  n  cos  a  +  6  sin  a 

Conversely, 

fi  =  i*  sin  a  +  3*  cos  a,  6  =  -1*  cos  a  +  j*  sin  a 

Differentiation  of  the  second  of  Eqs.  (1.6)  yields,  with  the  help  of 
Eq.  (1.4), 

-  t)sin  a  +  j*(|j  -  x)cos  a  +  k*(K*  cos  a  +  k  sin  a) 
Also,  Eqs.  (1.2)  and  (1.6)  yield 

3b  1  A  A 

*  -  |-(i*  sin  a  +  3*  cos  a) 

Equations  (1.7)  and  (1.8)  yield 


dot  1 
3s  E 


tan  a  =  - 


Also,  since  k*  =  t,  Eqs.  (1.2)  and  (1,4)  yield 


1-1  .  A 

3t  n  .  -> . 

31s  R  "  1  K  '  >** 


(1.5) 


(1.6) 


(1.7) 


(1-8) 


(1-9) 


(1.10) 


Equations  (1.6)  and  (1.11)  yield 

cos  a  ,  sin  a 
K  =  -  R  >  *'  *  -R- 


(1.11) 


(1.12) 


At  a  point  where  a  concentrated  couple  is  introduced  into  the  tube, 
classical  oeam  theory  indicates  that  the  curvature  1/R  is  discontinuous. 
Discontinuities  also  may  appear  in  1/E  ,  A,  and  6.  The  preceding  theory  is 
limited  to  points  at  which  the  pertinent  geometric  quantities  are  continuous. 


Figure  3.  Vectors  ?*,  j*,  n,  6  in  a  Normal  Plane  of  Curve  C* 

1.4  KINEMATICS  OF  THE  TUBE 

The  curve  C*  representing  the  deflected  axis  of  the  tube  at  time  t  is 
defined  by  the  vector  equation  r  *  r(s,t),  in  which  s  is  arc  length  on  the 
axis  of  the  tube  and  r  is  a  radius  vector  from  a  fixed  origin  to  the  point, 
s  on  the  curve.  The  vector  3r/3s  is  the  unit  tangent  t  of  the  axis  of  the 
tube.  The  vector  3r/3t  is  the  velocity  of  the  center  of  the  cross  section 
of  the  tube  at  point  s. 

The  triad  of  unit  vectors  (i*,  j* ,  k*)  may  be  conceived  to  be  glued  to 
a  cross  section  of  the  tube,  with  its  origin  at  the  center  of  the  cross 
section.  As  the  tube  deflects  and  twists,  that  cross  section  and  the 
attached  triad  (i*,  j*,  k*)  rotate  with  angular  velocity  ii)(s,t).  The  vector 

—  AAA 

w  is  resolved  into  components  ,  ui„,  in  the  directions  (i*,  j*,  k*) ; 
i .e.  , 

U)  =  1*0^  +  +  k*u>3  (1.13) 

AAA 

Since  (i*,  j*,  k*)  are  unit,  vectors  attached  to  the  tube,  the  following 
kinematic  relations  exist: 


/)'1  ^  •-  A  A  A 

-  to  x  j*  -•  k*;^  -  i*u3 

A 

a  to  X  Ic*  =  i*U>2  -  j*0)j  (1.14) 


Equation  (1.14)  is  similar  to  Eq.  (1.4),  but  the  physical  interpretation  is 
different.  Since  k*  -  t,  Eq.  (1.14)  yields 


O) 


1 


A 


(1.15) 


The  components  of  to  in  the  directions  fi  and  6  are  oo  =  w  •  ft  and 
a)^  =  w  •  6.  By  Eqs.  (1.6)  and  (1.13), 

w  =  w.  sin  a  +  w.,  cos  a  (1.16) 

n  1  2 


By  Eqs.  (1.6),  (1.14),  and  (1.16),  w  =  -b *  8t/9t.  Similarly,  to^  is  derived. 
Accordingly, 


to 


n 


(1.17) 


Subsequently,  <n_  is  designated  as  to  ,  in  which  the  subscript  "a" 

<1 

indicates  the  axial  component.  The  subscript  "t"  is  reserved  to  indicate 
the  partial  derivative  with  respect  to  time. 


1.5  VELOCITY  OF  THE  PROJECTILE 

Balloting  is  not  considered.  At  any  instant,  the  axis  of  the  projectile 
is  assumed  to  be  tangent  to  the  axis  of  the  deflected  tube  at  the  point 
where  the  centroid  of  the  projectile  lies. 

The  location  of  the  centroid  of  the  projectile  at  time  t  is  specified 
by  s  =  £(t) ,  where  £(t)  is  regarded  as  a  given  function.  The  absolute 
trajectory  of  the  centroid  of  the  projectile  accordingly  is  represented  by 

r  =  r[£(t) ,  t]  (1 -18) 

The  absolute  velocity  of  the  centroid  of  the  projectile  is  V  ■-  dr/dt,  where 
d/dt  denotes  the  total  derivative.  By  the  chain  rule  of  partial 
differentiation, 


E 


v  9r  df;  .  9r 

dt  "  3s  dt  3t 

Since  3r/3s  =  t,  this  yields 


V  = 


•  A 

St 


+ 


3r 

3t 


(1.19) 


Equation  ( 1 ^ 19)  signifies  that  the  absolute  velocity  of  the  centroid  of  the 
projectile  is  the  vector  sum  of  the  velocity  relative  to  the  contiguous  part 
of  the  tube  and  the  velocity  of  the  center  of  the  cross  section  of  the  tube 
at  which  the  centroid  of  the  projectile  lies.  This  conclusion  could  have 
been  anticipated  from  general  kinematical  theory. 

The  distinction  between  the  total  derivative  and  the  partial  derivative 
applies  to  any  function  of  s  and  t;  i.e., 


d(-) 

dt 


a(-)  .  h  aj^i 
at  ^  ‘  3s 


(1.20) 


where  (-)  denotes  any  function  of  s  and  t. 

In  view  of  Eq.  (1.19),  the  square  of  the  speed  of  the  projectile  is 


V2  =  l2  +  2£t  •  +  (|~)2 


(1.21) 


By  Eq.  (1.19),  the  component  of  velocity  of  the  projectile  tangent  to  the 
deflected  axis  of  the  tube  is 


V  =  V  •  t 
a 


K  + 


A 

t 


3r 

at 


(1.22) 


The  components  of  velocity  of  the  projectile  in  the  directions  of  the 
principal  normal  and  the  binormal  to  the  deflected  axis  of  the  tube  are 


V  •  fi  =  fi 


V 


u  =  b 


(1.23) 


It 


By  definition , 


~r2  -t  -t  .2 

A  -  A  •  A  =  A 


in  whieh  A  is  any  veotov. 
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Equation  (1.23)  signifies  that  any  component  of  velocity  of  the  projectile 
normal  to  the  axis  of  the  tube  is  the  same  as  the  corresponding  normal 
component  of  velocity  of  the  contiguous  tube. 

1.6  ACCELERATION  OF  THE  PROJECTILE 

The  acceleration  of  the  centroid  of  the  projectile  is 

~  _  dV  ?  9V  ,  3V 
a  '  dt  "  9s  3t 

By  Eqs.  (1.2)  and  (1.19), 


av  a  ft  at 

3s  “  K  R  +  3t 


SH 


av  ,  f{ .  i  M .  i*- 

3t  5t  «  9t  ,t2 


>«? 


Consequently, 


or,  since  3t/3t  =  w  x  t  (see  Eq.  (1.14)), 


•  2  2 — 
_  i”  ^  ••a.  v* 

a  =  n  ~  +  2£w  x  t  +  +  — s* 

R  3t2 


s=C 


(1.24) 


(1.25) 


The  terms  on  the  right  side  of  Eq.  (1.25)  can  be  identified  as  follows: 

(a)  The  centripetal  acceleration  of  the  centroid  of  the  projectile  relative 
to  the  momentary  form  of  the  axis  of  the  tube,  (b)  The  Coriolis  acceleration 
of  the  centroid  of  the  projectile.  (c)  The  tangential  acceleration  of  the 
centroid  of  the  projectile  relative  to  the  adjacent  part  of  the  tube. 

(d)  The  absolute  acceleration  of  the  center  of  the  cross  section  of  the 
tube  at  which  the  centroid  of  the  projectile  lies.  This  decomposition  of 
the  acceleration  could  have  been  anticipated  by  general  kinematical  theory. 

By  Eq.  (1.25),  the  components  a&,  a^,  a.D  of  the  absolute  acceleration 
a  in  the  directions  of  the  tangent,  the  principal  normal,  and  the  binormal 
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of  the  deflected  axis  of  the  tube  are 


A  —  **  A 

a=t*a=£+t 

a 


s=C 


•2  2— 
a  «  fi  •  a  =  —-  +  2^0)  ♦  6  +  n  • 
n  R  at2 


s-5 


O _ 


S 


6  *  a  =  -2£w  •  fi  +  S  •  ^-5* 

at 


(1.26) 


since 

n  •wxt  =  w*txn-w*6 

and 

A  —  ^  - .A  /\  . —  ^ 

b  *  in  x  t  =  u  •  t  x  b  =  -m  •  n 

1.7  ANGULAR  VELOCITY  OF  THE  PROJECTILE 

In  the  time  interval  dt,  the  cross  section  of  the  tube  at  which  the 
centroid  of  the  projectile  lies  undergoes  the  angular  displacement  to  dt . 

In  the  same  time  interval,  the  projectile  advances  the  distance  £  dt 
relative  to  the  tube.  The  angular  displacement  of  the  projectile  due  to  the 
latter  displacement  is  w£dt.  The  spin  of  the  projectile  relative  to  the 
tube  is  wt,  where  w  is  the  magnitude  of  the  spin.  (w  is  not  the  magnitude  of 
the  vector  <n.)  Consequently,  the  relative  angular  displacement  of  the 
projectile  during  dt  because  of  the  spin  is  wtdt.  The  absolute  angular 
displacement  of  the  projectile  during  dt  is  the  vector  sum  of  these 
components,  namely, 

wdt  +  £w  dt  +  wt  dt 

Therefore,  the  absolute  angular  velocity  of  the  projectile  is 

ft  =  w  +  £w  +  wt  (1*27) 


Equations  (1.3),  (1.13),  and  (1.27)  yield 


fi  a  W'-Dj  +  «£)  +  J*(W2  ♦  K'£)  +  k*(W3  +  *r£  +  U)) 


Now  to^  and  u>2  can  be  eliminated  from  Eq.  (1.28)  by  Eq.  (1.15) 
be  eliminated  by  F,q.  (1.12).  Thus,  Eq.  (1.28)  yields 


a  = 


* 


cos  a  +  J* 


sin  a) 


+  k*  (to  +  Wj  +  £t) 


By  the  vector-triple-product  theorem  (Ref.  1,  Appendix  C) : 


A  A  P)  t*  A  A  P)t*  A  A 

-t*j*  •  ||  *  jn*  •  ||  •  -  §  x  a-  x  j.) 


or,  since  i*  x  j*  =  k*  =  t, 


Bt  Bt  o  Bt 

^  +  l*i*  *  37  =  1  X  Bt 


Also,  by  Eq.  (1.6),  -l*  cos  a  +  J*  sin  a  =  6.  Consequently, 
yields 


— .  a  % 1  A  F  a  * 

Cl  =  t  x  +  b  ^  +  t(in  +  wa  +  £x) 


The  components  of  $1  in  the  fi,  6,  and  t  directions  are 


=  sT  ♦  h  ,  *»  a  •  6  ,  =  si  •  t 


Since  the  terms  in  the  scalar  triple  product  can  be  permuted 

*  -  a£  Bt  .  f  e  Bt 

n  ‘  1  x  at  -  H  '  n  x  1  -  ‘  n 


it  ft-  3; 


Equations  (1.30),  (1.31),  and  (1.32)  yield 


n„  -  -6  '  f!  ■  «b  ”  “  •  ll  *  I  ■  «a  =  “  f  “a  +  & 


(1.28) 
;  k  and  k'  can 

(1.29) 


Eq,  (1.29) 

(1.30) 

(1.31) 
cyclically, 

(1.32) 

(1.33) 
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liquations  (1.17)  and  (1.33)  yield 


n 


n 


(0 


% 


fla  ■  to  *  o>a  + 


(1.34) 


1.8  KINETIC  ENERGY  OP  A  BALANCED  PROJECTILE 

In  the  theory  of  kinematics,  any  reference  frame  is  admissible.  For 
kinetics,  however,  a  Galilean  reference  frame  must  be  introduced. 
Consequently,  the  vector  r  is  now  considered  to  be  specified  with  respect 
to  a  Galilean  reference  frame;  e.g.,  the  earth. 

The  kinetic  energy  of  translation  of  the  projectile  is 

1  —2 
T.  =  -  m  V“ 

V  L 


where  m  is  the  mass  of  the  projectile.  Consequently,  by  Eq.  (1.21), 


^  1  rt2  Dr  ,Di\2,  1  ,di\2 

Tt  -  2  +  ‘  at  +  ^Dt  ^  =  2  m  dt 


(1.35) 


The  kinetic  energy  of  rotation  of  the  projectile  is 

Tv  ■  T  b®,.2  *  O  *  I  l3  «a2 


(1.36) 


where  ij  is  the  moment  of  inertia  of  the  projectile  about  a  transverse  axis 
through  its  center  of  mass  and  is  the  moment  of  inertia  of  the  projectile 
about  its  longitudinal  axis.  The  total  kinetic  energy  of  the  projectile  is 
T  =  Tj.  +  Ty.  Consequently,  by  Eqs.  (1.33),  (1.35),  and  (1.36) 


Tl> 


♦  A 

2Ct 


Dr 

Dt 


-  + 


(~)2] 
^Dr  J 


\  ld(G 


Dt? 

Dt' 


(n 


Dt.  2 
Dtj 


+  *> 
•  •> 


+  2  V"  +  wa  + 


)' 


(1.37) 


Since  h  •  Dt/Dt,  $  ♦  Dt/Dt,  and  t  ♦  Dt/Dt  are  three  orthogonal 

A 

components  of  Dt/Dt, 


„a  Dt.  2  -  Dt.  2  ,  -  Dt .2  ,Dt  ? 

(b  ’  DtJ  +  tu  *  (^r} 


Dt 


DtJ 


"Dt 
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A  -  A 


Furthermore,  since  t  •  t  =  1,  t  *  3t/3t  =  0.  Consequently.  liq.  (1.37) 
reduces  to 


Tp  *  7  *  &  ■  &  *  <!t>2i  ♦  I  V<H>2  * 2 1 *  •  I  *  S’ 


"  7  V"  *  “» * ,"'2 


(1.38) 


A  A 

Since  3t/Ds  =  n/R  and  1/R“  ■-  (3t/9s)  ,  Bq.  (1.38)  may  be  expressed  as 
follows : 


rr  „  1  r Sr  ,  ,3i\2.  1  .  ,3t  2 

Tp  =  j  m[4  -  25t  *  ^  "■  (w)  ]  +  -  i,(^r  +  C  — ) 


2  P3t  ^  3sJ 


+  \  V"  +  wa  "  ^T)2 


(1.39) 


or  more  concisely, 


tp  ■  7  ”(3F>2  *  7  b't*2  *  7  *  «U 


-•_,2 


(1.40) 


1.9  RELATION  BETWEEN  TWIST  AN!)  ANGULAR  VELOCITY 

Torsional  vibrations  of  a  perfectly  straight  tube  exhibit,  a  simple 
relationship  between  twist  and  angular  velocity.  During  a  time  interval  dt, 
the  cross  section  of  the  tube  at  point  s  undergoes  the  angular  displacement 
a)  dt.  During  the  same  time  interval,  the  cross  section  at  point  s  +  ds 
undergoes  the  angular  displacement 

3w 

(w  +  ™-  ds)dt 
■  a  os 

Consequently,  the  increment  of  twist  at  point  s  during  dt  is  (3w  /3s) dt. 

A 

Also,  si, ice  the  total  twist  at  point  s  is  x  -  Tp,  the  increment  of  twist  at 
point  s  during  dt  is  (3x/3t)dt.  Therefore,  for  a  straight  tube  that 
executes  torsional  motion, 


Equation  (a)  might  be  adopted  as  an  approximation  for  a  slightly  curved 
tube,  but  an  inconsistency  arises.  It  is  illustrated  by  a  perfectly  rigid 
curved  tube.  In  this  case,  t  =  Tq,  since  there  is  no  deformational  twist  of 
a  rigid  tube.  Accordingly,  Eq.  (a)  yields  8to  /8s  =  0.  Therefore,  to  does 

ft  cl 

not  depend  on  s.  Furthermore,  since  all  particles  of  a  rigid  body  have  the 
same  angular  velocity  at  any  instant,  the  vector  angular  velocity  w  of  a 
rigid  tube  is  independent  of  s.  The  axial  component  of  angular  velocity  of 
a  cross  section  is  to  -  to  •  t.  For  a  curved  tube,  t  obviously  depends  on  s. 
Thus,  we  arrive  at  the  contradictory  conclusions  that  to  depends  or  s  and 

cl 

to  does  not  depend  on  s.  Consequently,  in  general,  Eq.  (a)  must  be  rejected, 
a 

although  it  is  correct  in  certain  special  cases. 

Since  a  theory  of  flexible  tubes  must  be  consistent  with  the  theory  of 
rigid-body  displacements,  the  nature  of  the  function  to  (s,t)  for  a  rigid 

cl 

curved  tube  is  pertinent.  In  general,  to  =  to  •  t,  and,  for  a  rigid  tube, 
to  =  w(t) .  Therefore,  for  rigid  tubes. 


In  view  of  Eqs.  (1.1)  and  (1.2),  this  yields 


*  1  —  ^  A 
n  =  jj-i.o»bxt 


Since  the  vectors  in  the  scalar  triple  product  may  be  permuted  cyclically, 
this  yields 


vu.i  1  * 

a  1  C  ”  2 

3T'  -  Rb  '  "x  ' 


Since  to  x  t  -  9t/8t,  this  yields  the  following  equation  for  rigid  tubes : 


a  l  f  3t 

sT  *  if  b  •  at  * 0 


The  general  theory  of  flexible  curved  tubes  is  now  considered.  An 
infinitesimal  segment  of  such  a  tube,  as  viewed  along  the  binormal  of  the 
axis  is  shown  in  Figure  4.  Infinitesimal  angular  displacements  are 


resolved  along  the  horizontal  line  L.  To  first-degree  quantities,  the 
component  on  line  L  of  the  rotation  of  the  right-hand  cross  section  relative 
to  the  left-hand  cross  section  is 

-~-^.ds  dt  +  to,  d@  dt 
9s  1 

Accordingly,  the  increment  of  twist  at  section  s  during  dt  is 


/a 

(9s  +  W1  9s^dt 


Since  this  is  equal  to  9t/9t  dt, 


9t  3a>a  36 

9t  9s”  W1  3s 


Figure  4.  Resolution  of  Angular  Velocities 


With  the  sense  of  indicated  in  Figure  4,  =  -o)  •  n.  Also,  30/ 3s  =  1/R. 

Therefore , 


3t  _  __a  a)  •  n 

3t  3s  R 

Since  n  =  f>  x  t ,  this  yields 


3t  ^Wa  1  —  £  £  1  £  —  •*> 

w*  w  -  r“  • b  x  * =  sr 4  r  b  •  “ x  1 


Since  u)  x  t  =  3t/3t,  this  yields 


3r  ^_a  l_fi  3t 

at  “  as  ‘  +  r  b  *  at 


(1.41) 


For  torsional  motion  of  a  straight  tube,  Eq.  (1.41)  yields  Eq.  (a) 
since  1/R  =  0.  For  motion  of  a  rigid  tube,  it  yields  Eq .  (b) ,  since 
3r/3t  =  0.  If  the  axis  of  the  tube  bends  in  a  plane  F,  vector  b  is 
perpendicular  to  F,  and  vector  3t/3t  lies  in  plane  F.  Consequently, 

S  -  3t/3t  -  0.  Accordingly,  if  curves  C  and  C*  are  constrained  to  lie  in  a 
plane,  Eq.  (1.41)  reduces  to  Eq.  (a).  Usually,  in  studies  of  plane  motion 
of  gun  tubes,  x  and  are  taken  to  be  zero. 

a 

Equation  (1.41)  can  be  expressed  more  simply,  since 


1  c  ot  1  o  —  £  1  —  *  c 

„6.^=f6.U,Xt=^U3-tx6  = 


i  -  -  at 

_  w  .  n  -  _w  . 


Also, 


:(u)  *  t)  -  03 


Therefore,  Eq.  (1.41)  reduces  to 


(1.42) 


For  a  rigid  tube,  3i/3t  -  0  and  u)  =  w(t)  ,  so  Eq.  (1.42)  is  satisfied. 
Equation  (1.42)  means  that  3t/3t  is  equal  to  the  tangential  component  of 
3o)/  3s . 


In  view  of  Eq.  (1.17),  Eq.  (1.41)  may  be  expressed  as  follows: 

9(0  to 
dx  an 


at  9s  R 


(1.43) 


The  general  solution  of  Eq.  (1.43)  is  easily  derived.  There  is  a  function 
i|i(s,t)  such  that  t  =  9)^/9s.  The  function  tp  contains  an  arbitrary  additive 
function  of  t.  Also,  there  is  a  function  X(s,t)  such  that  w  =  X  +  3\J>/9t. 

9 

Equation  (1.43)  accordingly  yields 


(0 

3a  _  n 

as  ‘  r 


Consequently, 


'S  0) 

X  =  ds  + 

Jo  R 


The  arbitrary  additive  function  of  t  in  the  function  ip  may  be  chosen  to 
cancel  q(t)  .  Therefore,  there  is  a  function  i|i(s,t)  such  that 


9iji  9il)  fs  wn  . 

9s  ’  a  3t  L  R 


(1.44) 


The  function  ip (s , t )  represents  the  angular  displacement  of  a  cross  section 
of  the  tube  in  its  plane. 

Equations  (1.34)  and  (1.44)  yield 


£2  =  co  +  u)  +  £t  =  u>  + 

a  a 


w  , .  s=C 


(1.45) 


Equation  (1.45)  may  be  substituted  into  the  kinetic  energy  expression  for 
the  projectile  (Eq.  (1.39)  or  (1.40)). 


1.10  VIRTUAL  WORK  ASSOCIATED  WITH  THE  PROJECTILE 

In  order  to  apply  Hamilton's  principle  or  Lagrange's  equations  to 
problems  of  gun  dynamics,  we  require  the  expression  for  the  virtual  work  of 
all  the  non-inertial  forces  that  operate.  It  is  a  linear  expression  in  the 


infinitesimal  virtual  displacements.  Only  the  part  of  the  virtual  work  that 
explicitly  involves  the  projectile  is  considered  in  this  section. 

Conceptually,  the  system  receives  an  infinitesimal  virtual  displacement 
that  generally  does  not  coincide  with  the  true  course  of  the  motion.  The 
real  motion  of  the  system  is  imagined  to  be  stopped  while  the  virtual 
displacement  is  performed.  The  actual  forces  in  the  system  are  imagined  to 
persist  while  the  virtual  displacement  is  executed.  In  the  present  case, 
the  vector  r(s,t),  defining  the  curve  C*,  receives  a  virtual  increment  6r(s). 
This  transforms  curve  C*  into  another  curve  C**,  Figure  5.  Since  the  center 
line  of  the  tube  is  considered  to  be  inextensional ,  the  variation  <5r  must 
conform  with  this  constraint.  Consequently,  the  mapping  C*  ■+  C**  must  be 
performable  by  inextensional  bending  and  twisting  of  curve  C*.  In  addition 
to  the  virtual  displacement  5r,  the  coordinate  £  of  the  projectile  receives 
a  virtual  increment  S£,  and  it  also  receives  a  virtual  angular  displacement 
about  its  longitudinal  axis. 


Figure  5.  Virtual  Displacements 
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The  work  that  a  force  F  perforins  on  a  particle  during  a  time  interval 
(tQ,  t ^ )  is  defined  by 


W 


V  dt 


where  V  is  the  velocity  of  the  particle.  Since  V  depends  on  the  choice  of 
the  reference  frame,  so  does  W;  i.e.,  work  is  a  relative  quantity.  This 
conclusion  is  consistent  with  the  fact  that  the  work  of  all  the  forces  that 
act  on  a  system  equals  the  increase  of  kinetic  energy  of  the  system,  since 
kinetic  energy  also  is  a  relative  quantity,  inasmuch  as  it  depends  on  the 
velocities  of  the  particles.  In  this  analysis,  work  is  calculated  with 
respect  to  the  Galilean  reference  frame  to  which  the  vector  r  is  referred. 

The  contact  forces  that  the  projectile  exerts  on  the  bore  are  reacted 
by  equal  and  opposite  contact  forces  that  the  bore  exerts  on  the  projectile. 
Consequently,  the  normal  components  of  all  these  forces  perform  no  net  work 
on  the  system.  Accordingly,  the  gyroscopic  couple  of  the  projectile  and  its 
reaction  perform  no  net  work.  An  analogy  is  a  person  who  lifts  an  object. 
The  person  performs  work  on  the  object,  but  gravity  performs  an  equal  amount 
of  negative  work.  Together,  the  lifter  and  gravity  perform  no  net  woi-k. 

The  projectile  receives  the  virtual  displacement  <$£  relative  to  the 
tube,  and  the  contiguous  tube  receives  the  axial  virtual  displacement 

a  _ 

t  •  Sr.  Consequently,  the  absolute  axial  component  of  the  virtual 
displacement  of  the  projectile  is  +  t  •  <Sr.  The  driving  force  on  the 

base  of  the  projectile  is  P^A,  where  is  the  pressure  of  the  gas  and  A  is 
the  cross  sectional  area  of  the  bore.  There  is  a  resisting  force  F  which 
results  from  axial  friction  and  engraving  of  the  rifling.  Also,  there  is  a 
resisting  force  P2A  resulting  from  pressure  of  the  air  ahead  of  the 
projectile.  Consequently,  the  virtual  work  of  forces  that  act  on  the 
projectile  is 

(P:A  -  P2A  -  F)(«£  +  t  •  6r) 

/N  — 

The  virtual  work  of  friction  and  engraving  on  the  tube  is  Ft  *  or.  The  net 
virtual  work  6W1 ,  resulting  from  axial  movement  of  the  projectile,  is  the 
sum  of  these  expressions.  Consequently, 
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I  . 


n«»  V-#  ’ sTf 7 •t'*Wi,~J >*■  ■  U"‘ • 


6W1  =  [(P1  -  P2)A  -  F]<$5  +  (Pj  -  P2)A  t  •  6r 

in  which  Sr  is  evaluated  at  the  point  s  -  £. 

Also,  there  is  a  contribution  6W^  to  the  virtual  work  from  the  rifling 
torque.  The  virtual  angular  displacement  of  the  projectile  relative  to  the 
tube  is  5x,  where  x  =  w.  The  angular  displacement  of  the  cross  section  of 
the  tube  at  point  s  =  5  is  denoted  by  \l>.  The  absolute  virtual  angular 
displacement  of  the  projectile  is  <$x  *  <$’/'•  The  rifling  torque  is  denoted  by 
M , .  The  virtual  work  performed  on  the  projectile  by  the  rifling  torque  is 
MpC^X  +  S40  •  The  virtual  work  performed  on  the  tube  by  the  rifling  torque 
is  -Mf  6ip.  Consequently, 

6W2  =  Mf(6X  +  6<p)  -  Mr  6<p  =  Mr  SX 

Aside  from  effects  of  gravity,  the  virtual  work  explicitly  related  to 
the  projectile  is  SlV  =  SWj  +  S IV^.  Hence, 

6W  =  [ (PL  -  P2)A  -  F]$£  *  (Pj  -  P2)A  t  *  67  +  Mr  6X  (1.46) 

in  which  relevant  functions  are  evaluated  at  the  point  s  =  £.  Additional 
contributions  to  the  total  virtual  work  of  the  system,  coming  from  the 
action  of  gas  pressure  on  the  breech,  effects  of  gravity,  strain  energy  of 
the  tube,  and  effects  of  the  supporting  structure,  are  not  considered  here. 

We  adopt  the  viewpoint  that  ?(t)  and  u> ( t )  are  given  functions.  Then 
S£  =  <5x  =  0,  and  Eq.  (1.46)  is  simplified  accordingly.  Also,  Sr  is 
restricted  by  the  condition  of  inextensionality  of  the  tube.  Since 
t  •  t  =  1,  t  •  6t  =  0.  Consequently,  since  t  =  3r/9s, 

t  •  6rs  =  0  (1 .47) 

Equation  (1.47)  expresses  the  constraint  on  Sr. 

If  the  tube  is  initially  straight,  t  is  approximately  a  constant 
vector,  since  the  deflections  are  small.  With  this  approximation, 

Eq.  (1.47)  yields 

^  ^  _ 

-5— (t  •  Sr)  =  0  or  t  •  Sr  =  constant 
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If  r  is  given  at  the  breech  (s  =  0) ,  Sr  *  0  at  the  breech.  Then,  since 
—  /\  — 

t  ♦  5r  =  constant,  t  *  or  =  0  everywhere.  Accordingly,  for  a  gun  with  a 
straight  tube,  Eq.  (1.46)  yields  <5W  =  0.  The  condition  t  •  5r  =  0  means 
that  6r  must  be  perpendicular  to  the  axis  of  the  bore. 


1.11  KINEMATIC  .AND  GEOMETRIC  RELATIONS  IN  SCALAR  NOTATION 

Adaptation  of  the  foregoing  theory  to  digital  computer  programming 
requires  that  the  equations  be  expressed  in  scalar  form.  Rectangular 
coordinates  (x,  y,  2)  with  corresponding  unit  vectors  (1,  j,  tc)  are 
attached  to  a  Galilean  reference  frame,  but,  insofar  as  kinematics  is 
concerned,  the  reference  frame  is  arbitrary.  The  orientation  of  these  axes 
with  respect  to  the  gun  is  not  of  immediate  concern.  The  deflected  axis  C* 
of  the  tube  at  time  t  is  defined  by  the  equations,  x  =  x(s,t),  y  =  y(s,t), 
z  =  z(s,t).  Since  s  is  defined  to  be  arc  length  on  curve  C* , 

2  2  2  , 

Xs  +  ys  +  Zs  =  1  (1-48) 


in  which  the  subscript  denotes  the  partial  derivative;  e.g.,  x  =  3x/3s. 

s 

The  radius  vector  from  the  origin  of  the  (x,  y,  z)  coordinates  to  a  point 
on  curve  C*  is 

r=ix+jy+kz  (1.49) 

The  unit  tangent  of  curve  C*  is 


/'N 

t 


3r 

3s 


^  A 

ix  +  1  y  +  kz 
s  J  '  s  s 


(1.50) 


By  Eqs .  (1.2)  and  (1.50), 
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II 
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IX  + 


ss 


(1.51) 


Therefore, 


By  Eqs.  (1.1)  and  (1.51),  S/R  may  be  expressed  in  determinant  notation,  as 
follows: 
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b  _  * 
R  =  1 


By  Eq.  (1.50), 

9t  ?  .  ❖ 


dt  ^  ^ 

37  =  1  Xst  +  J  yst  +  *  zst 


(1.53) 


(1.54) 


By  Eqs.  (1.17),  (1.51),  (1.53),  and  (1.54), 


—  =  x  x  +  y  y  +  z  z  _ 
R  ss  st  yss  'st  ss  St 


(1.55) 


U>n  xst  ^st  2st 

T=  xs  ys  zs 

xss  yss  zss 


(1.56) 


Accordingly,  and  oj^  are  determined,  if  the  functions  x(s,t),  y(s,t), 
z(s,t)  are  known.  The  component  w  is  not  determined  solely  by  these 
functions,  since  it  depends  on  the  twist  of  the  tube,  as  is  indicated  by 
Eq.  (1.43). 

The  tortuosity  of  curve  C*  is  determined  most  readily  by  the  second  of 
Eqs.  (1.2).  It  yields 


1  -  3fi 

r  •  b  •  sf 


(1.57) 


By  the  first  of  Eqs.  (1.2), 


3  t  _3A  1  W  .  a  All 

.  2  "  3sl‘RJ  "  R  3s  n  3s 
3s 


f-  ? 

ix  +  ]  y  +  k  z 
sss  sss  sss 


(1.58) 
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Therefore, 


£  .  i_A  „  R-2g  .  M  .  _i_ 

R  SsW  0  9s  r2j. 


With  Eqs.  (1.53)  and  (1.58),  this  yields 
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(1.59) 
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It  is  possible  to  eliminate  R  from  Eq.  (1.59)  by  means  of  Eq.  (1.52).  Thus, 
l/l  is  expressed  as  a  rational  function  of  first,  second,  and  third 
derivatives  of  x,  y,  and  z  with  respect  to  s. 


1.12  VELOCITY  AND  ACCELERATION  OF  THE  PROJECTILE  IN  SCALAR  NOTATION 

The  velocity  of  the  centroid  of  the  projectile  is  given  by  Eq.  (1.19) 
Consequently,  by  Eqs.  (1.49)  and  (1.50), 


=  +  1  xs)  +  J(yt  +  k  ys)  +  k(zt  +  i  zs) 


S  =  5 


or 


V 


?  dx  »  dy  c 
1  7r  +  1  T r  +  K 
dt  J  dt 


dz 

dt 


s-£ 


where  d/dt  denotes  the  substantial  derivative.  Hence, 


(1.60) 


=  A2  + 

(dt7  + 


$>2 


(-V 

'■dt' 


s=S 


(1.61) 


The  axial  component  of  velocity  of  the  projectile  is 


V 

a 


=  V 


■■  dx 

t  -  X  y-  +  y 
c  dt  7 


dy 

dt 


+  z 


dz 

dt 


s=C 


(1.62) 


The  components  of  V  on  the  principal  normal  and  the  binormal  of  the 
deflected  axis  C*  are  given  by  Eq.  (1.23).  Consequently,  by  Eqs.  (1.49), 
(1.51) ,  and  (1.53)  , 


36 


(1.63) 


x  x  +  y  y.  +  z  z. 
SS  t  '  ss  't  ss  t 
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(1.64) 


The  acceleration  of  the  centroid  of  the  projectile  is  given  by 
Eq.  (1.24).  Consequently,  by  means  of  Eqs.  (1.49),  (1.50),  (1.51),  and 
(1-54), 


a  =  i(£2  x  +  2  f ‘x.+£x  +  x.  .) 

v  ss  st  s  tt' 


.  Set2  yss  ♦  2  t  yst  •  5  y,  *  ytt) 


+  £(£2  z  +  2  5  z  .  +  E  z  +  z> 

v  SS  ^  St  ^  S  tt' 


S=C 


(1.65) 


Therefore,  by  Eq.  (1.26),  the  axial  component  of  acceleration  of  the 
centroid  of  the  projectile  is 


a  =  t  *  a  =  £  +  X  X  +  y  y  +  z  z.  . 
a  *  s  tt  's  'tt  s  tl 


(1.66) 


By  Eq.  (1.26),  the  component  of  a  on  the  principal  normal  vo  curve  C* 
is  determined  by 


A  <i 

a  *  n  _  ji 
R  "  "R 


•  2  •  vs.  2— 

Z *  ,  K  A  fi  3  r 

r2  2RWb  R  *  3,2 


s=^ 


Consequently,  in  view  of  Eq.  (1.51), 


a  i2  t 

R  „2  +  2  R  wb  +  xss  Xtt  +  ^ss  ^tt  +  Zss  Ztt 
K 


(1-67) 


The 


factor  w^/R  in  Eq.  (1.67)  can  be  eliminated  by  means  of  Eq.  (1.55). 


37 


Likewise,  the  component  of  a  in  the  direction  of  the  binormal  b  of  curve 
C*  is  determined  by  Eq.  (1.26).  In  view  of  Eq.  (1.53), 

ytt  *tt 

ys  zs  (1.68) 

y  z 

ss  ss 

The  factor  w^/R  in  Oq.  (1.68)  can  be  eliminated  by  means  of  Eq.  (1.56). 

Also,  R  can  be  eliminated  by  means  of  Eq.  (1.52). 


1.13  KINETIC  ENERGY  OP  THE  PROJECTILE  IN  SCALAR  NOTATION 

The  kinetic  energy  of  the  projectile  is  given  by  Eq.  (1.40). 
Consequently,  by  Eq.  (1.50), 

1  r,dx.  2  .dy.  2  .  .dz.2, 

V?',l<af)  ♦  (jr  *  <jf>  i 
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1  .  ■>  c=r 

+  y  i3Cw  +  w#  +  £t)_  "  ■* 


(1.69) 


1.14  APPROXIMATE  THEORY  FOR  INITIALLY  STRAIGHT  TUBES 

If  the  tube  is  initially  straight,  the  z-axis  is  conveniently  chosen 
to  coincide  with  the  undeflected  axis  of  the  bore.  Then  x(s,t)  and  y(s,t) 
are  deflection  components  of  the  tube.  By  Eq.  (1,48), 


=  Vl  - 


7  2  12  2 

(xs  +  v  =  1  ■  l(xs  +  ys > + 


(1.70) 


Since  the  deflections  of  the  barrel  of  a  gun  are  small,  it  is  reasonable  to 
approximate  Eq.  (1.70)  by  =  1.  Then, 


The  function  f(t)  represents  the  value  of  z  at  the  point  s  -  0.  This 
point  is  conveniently  chosen  to  lie  at  the  breech.  It  may  vary  with  time, 
if  the  breech  moves.  It  is  to  be  noted  that  F.q.  (1.71)  signifies  that 
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x  +  y,/  =  0.  Therefore,  the  accuracy  of  quadratic  expressions  in 

S3 

derivatives  of  x  and  y  requires  study.  In  view  of  Eq.  (1.71),  Eq.  (1.49) 
becomes 


r  =  .x  x  +  J  y  +  ft  s  +  fc  f(t) 


(1.72) 


Accordingly,  the  unit  tangent  of  curve  C*  la 


A  A  A 

t  =  1  xs  .  j  ys  . 


(1.75) 


By  Eq.  (1.51), 


ii  <•  * 

tr  =  ix  +  i  y 

R  ss  J  'ss 


(1.74) 


By  Eq.  (1.52), 


iT^ss  +yss 


(1.75) 


If  curve  C*  is  constrained  to  lie  in  the  yz  plane,  x  =  0  and  Eq.  (1.75) 
reduces  to  1/E  -  +_  ygs,  which  is  a  well-known  linear  approximation  in  the 
engineering  theory  of  beams.  However,  if  x$;.  and  yss  both  differ  from  zero, 
there  is  no  linear  approximation  of  Eq.  (1.75)  available  by  Taylor  series 
expansion. 

The  linear  approximation  of  6/R,  obtained  fi'om  Eqs.  (1.53)  and  (1.71), 


—  =  -ty  +  j  x 
R  • ss  '  ss 


(1.76) 


Here,  the  term  ^(x^  y  •  ys  xss)  has  been  discarded,  since  it.  presumably  is 
small  compared  to  the  linear  terms  in  Eq.  (1.76). 

In  view  of  Eqs.  (1.59)  and  (1.75),  the  tortuosity  of  curve  C*  is 
approximated  by 


XV  -  y  X 

1  ' ss  'sss  1  ss  sss 

Z  =  2“  2 

xss  +  ySS 


(1.77) 
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Equations  (1.55)  and  (1.50)  reduce  to 


- —  s  v  v  >  V  V 
K  SS  St  XSS  *  St 


R  °  Xst  ^  SS  ^  St  X3S 


(1.78) 


Introduction  of  the  foregoing  approximations  into  the  equations  for  the 
velocity  and  the  acceleration  of  the  projectile  (Art.  1.12)  is  routine. 
Equation  (1.69),  which  gives  the  kinetic  energy  of  the  projectile,  becomes 

Tp  ■  I  *  $>2  *  d  *  i-)2i 

(yst  *  <  >W>2' 


1  .  .  •  ..2  s=£ 

'•  2  X3  (w  +  %  +  ^ 


(1.79) 


By  F.qs .  (1.45)  and  (1.78) , 


ft  =  to  +  (o  +  £r 

tl  (l 


it  S~5  =  »  +  f  (x  .  yoc  -  y_  x  Ids  + 


st  7  ss  7 st  ss7'1"  dt 


(1.80) 


Consequently,  Eq.  (1.79)  yields 


T  “  \  >"[xt.2  +  yt2  +  2C(xs  xt  +  yg  yt)  +  (t  +  f)2] 


+  i-  i,  [x  2  +  y  2  «  2C(x  X  .  +  y  y,J  |2(x  2  +  y  2)] 


SS  St.  '  SS  St 


SS  SS 


I  Vw  +  *t  +  +  L  (xst.  yss  ■  yst  xss)tl5l 


2  s=C 


(1.81) 


2  2 

Since  the  approximation  ?.  =  s  is  used,  Eq.  (1.70)  implies  that  xg  +  yg  =  f 
so  this  expression  has  been  discarded  from  Eq.  (1.81).  Also,  because  of 
the  approximation  z  =  s,  terms  of  third  and  fourth  degree  in  Eq.  (1.81) 
have  little  credibility.  Consequently,  F.q.  (1.81)  is  simplified  by 


elimination  of  these  terms.  Thus,  Eq.  (1.81)  is  reduced  to  the  following 
quadratic  expression  in  x,  y,  i|i  and  their  derivatives: 


T 

P 


1  2  .  2 


2  m[xt  +  yt  +  2^(xs  Xt  +  ys  yt}  +  +  ^21 


s=$ 


1  .  ,  2  2 
+  i'st  * 


ss  st  'Ss  'stJ 


+  2  i3(w  +  +  ^s)2|S"5  +  ^4  (xst  yss  -  yst  xss)ds  n-82) 


j 


1 
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SECTION  2 

FORCES  AND  MOMENTS  ACTING  ON  A  BALANCED  PROJECTILE  IN  A  FLEXIBLE  TUBE 


2.1  INTRODUCTION 

Section  1  deals  primarily  with  kinematics  of  a  flexible  tube  containing 
an  accelerating  projectile.  In  this  section,  the  theory  is  extended  to 
provide  formulas  for  the  forces  and  moments  acting  on  the  projectile.  By 
specialization,  the  forces  and  moments  acting  on  a  projectile  in  a  moving 
rigid  tube  are  obtained. 

As  in  Section  1,  balloting  is  excluded.  The  weight  of  the  projectile 
is  disregarded,  but  it  would  merely  augment  the  forces  on  the  projectile  by 
the  term  mg,  where  g  is  the  vector  acceleration  of  gravity.  It  would  have 
no  effect  on  the  moment  vector,  except  indirectly,  through  its  dynamic 
effect  on  the  deflection  of  the  tube.  The  dynamic  response  of  the  tube  is 
not  treated  in  this  section. 

2.2  FORC"  ON  THE  PROJECTILE 

The  net  force  on  the  projectile  is  F  =  m  a.  Consequently,  if  the 
center  of  mass  of  the  projectile  coincides  with  the  centroid,  Eq.  (1.25) 
yields 


—  N  £  - -  A  #*  A  9  I* 

F  =  in(n  +  2£<o  x  t  +  E,t  +  — y) 

K  9t 

By  Eq.  (1.26),  the  n,  6,  t  components  of  F  are 


F 

a 


m(5 


+  t 


(2.1) 


F  =  m(i~  +  2£w  •  b  +  n  • 
n  K 

j.2— 

F,  =  m(-2?io  •  ft  +  $  •  5_T)  (2.2) 

0  9t 

wherein  relevant  functions  are  evaluated  at  point  s  =  E. 


The  axial  force  on  the  projectile  due  to  friction  and  engraving  is 
denoted  by  F.  Its  positive  sense  is  toward  the  breech.  The  pressure  on  the 
base  of  the  projectile  is  denoted  by  P^  and  the  resisting  pressure  of  air 
ahead  of  the  projectile  is  denoted  by  P^.  Accordingly, 

Fa  =  -F  ♦  (Pt  -  P2)A 

where  A  is  the  cross  sectional  area  of  the  bore.  Therefore, 

j.2— 

F  =  (Pj  -  P2)A  -  »(?  ♦  t  •  ^-|)  (2.3) 

3t 

2—  2 

The  factor  3  r/3t  is  the  acceleration  of  the  center  of  the  cross  section 
of  the  tube  at  which  the  center  of  mass  of  the  projectile  lies. 

2.3  MOMENT  OF  FORCES  ACTING  ON  THE  PROJECTILE 

Because  of  axial  symmetry  of  the  projectile,  two  of  its  principal 
moments  of  inertia  are  equal;  i.e.,  i^  =  i2-  Consequently,  the  components 
of  angular  momentum  of  the  projectile  with  respect  to  its  center  of  mass 
are 


H  =  i,  ft  , 
n  In 


Hb  "  ll  S  *  Ht  =  13  fia 


(2.4) 


Hence, 


H  =  i,(fi  ft  +  6  ft.  )  +  i,  t 
lv  n  b  3 


ft 


The  vector  ft  is  given  by  Eq.  (1.34) 
Since  3fi/3t  =  w  x  fi,  etc.. 


3n  —  a  a 

■a—  -  u>  x  fi  =  b  u  -  t  a), 
3t  a  b 


■  w  x  6  3  t  w  -  fioj 
3t  n  a 


3t  —  i  ^  C 
77--  -  in  x  t  =  n  u).  -  b  a) 

3t  b  n 


(2.5) 


(2.6) 


Equations  (2.6)  tire  analogous  to  Eq.  (1.14). 


Q 


The  moment  about  the  center  of  mass  of  the  projectile  of  the  forces 
that  act  on  the  projectile  is 


rr  dH  3H  i  JH 

M  =  dt  =  3t  +  C  a? 


The  pressure  force  (P^  -  P^jA  and  the  weight  exert  no  moment  about  the 

center  of  mass  of  the  projectile.  Consequently,  the  moment  M  results 

entirely  from  contact  forces  between  the  projectile  and  the  bore.  The 

rifling  torque  is  M  •  t  =  M  .  The  gyroscopic  couple  is  M  •  n  =  M  . 

st  n 

Equations  (2.5)  and  (2.7)  yield 

»  -  &(«  “n  *  6  V  *  b  jf(?  V 

The  derivatives  3fi/3s,  36/3s,  3t/3s  are  given  by  Eq.  (1.2).  Since  dft/dt  = 
3n/3t  +  E,  3ft/ 3s,  etc.,  Eqs.  (1.2)  and  (2.6)  yield 

dfi  f  f  i,6  t. 

If  ’  b  “a  -  *  “b  *  -  W> 


/V.  *  II 

ft  <*»a  -  C  J 


b  w  +  £,  ~ 
n  ^  R 


It  is  convenient  to  introduce  the  notation. 


(2.7) 


(2.8) 


(2.9) 


«  +  i  =  A 

a  Z 


Equations  (1.34),  (2.9),  and  (2.10)  yield 


(2.10) 


£  -  bX  -  t  SL  ,  £ 
dt  b  dt 


ta  -fix,  =  ft  a  -  6  ft 

n  dt  b  n 


(2.11) 


Consequently, 


~(n  fi  ) 
dt v  ir 


d  ft 

fi  •  +  ft  (6a  -  t  a. ) 
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(2.12) 


d  ft 

--(6  ft.)  =/b  .  ... 

dt  b'  dt  n 


+  ft,  (t  ft  -  fU) 


d  ft 

-—(£  ft  )  =  t  +  ft  (fi  ft.  -  6  ft  ) 

dt v  a'  dt  b  n' 


Equations  (2.8)  and  (2.12)  yield 

„  =  yfi  ♦  6  -  XCfi  fl„  -  6  «„)]  .  ljt«  ^  *  9aC*  S  -  6  V) 

(2.13) 


Hence, 


d  ft 


M  =  i.  — j—  +  ft,  (i_  ft  “  i-i^) 
n  1  dt  bv  3  a  1  ’ 


\  - 


=  i 


.  dS> 


-  n„(i,  ft-  -  M) 


1  dt  aV*3  a  1 


d  ft 


M  *  i 


a  3  dt 


(2.14) 


Equations  (2.14)  resemble  the  Euler  equations  (with  i^  s  i2) ,  but  they 
are  not  identical  to  them  unless,  by  chance,  ft  -  The  explanation  for 
the  difference  lies  in  the  physical  reference  frames  to  which  the  equations 
refer.  Euler's  axes  (1,  2,  3)  are  the  principal  axes  of  inertia  of  the 
projectile.  They  are  imbedded  in  the  projectile.  For  example,  Euler’s  ft^ 
is  the  component  of  ft  on  a  transverse  axis  that  is  imbedded  in  the 
projectile.  On  the  other  hand,  ft^  is  the  component  of  ft  on  the  principal 
normal  of  the  trajectory  of  the  center  of  mass  of  the  projectile.  At  a 
given  instant,  these  two  axes  may  coincide,  so  that  ft^  =  ft^,  but  the  time 
rate  of  change  of  ft.^ ,  denoted  by  d  ft^dt,  is  generally  not  the  same  as  the 
time  rate  of  change  of  ft^.  The  time  rate  of  change  of  the  orthogonal 
projection  of  vector  ft  on  a  moving  axis  clearly  depends  on  the  motion  of 

a 

that  axis,  and  the  vectors  fi,  b,  do  not  have  the  same  motions  as  the  lateral 
principal  axes  of  inertia  of  the  projectile. 
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2.4  MOMENTS  ACTING  ON  A  PROJECTILE  IN  A  MOVING  RIGID  TUBE  OF  ANY  FORM 
If  the  tube  is  rigid,  its  angular  velocity  is  io(t) ;  i.e.,  w  is 
independent  of  s.  Also,  the  deformational  twist  x  is  zero.  Furthermore , 
R  =  R(s)  and  E  =  E(s).  Equations  (1.2),  (1.34),  and  (2.10)  apply  without 
alteration  (except  that  x  =  0) .  Since  w  =  to(t)  and  id  =  w  •  n,  etc., 


3w  (o.  a)  3to, 

_ n  —  3h  b  a  _l 

3s  ‘  W  *  3s  =  E  "  R  *  3s 


(2.15) 


Therefore, 


do)  3  w  0).  w  do), 

n  n  t,  b  a.  b 

dt  "  3t  ^ E  '  R  J  '  dt 


30),  •  do) 

_ b  £  _ a 

3t  "  I  wn  *  dt 


K  ^  i 

sr +  R  “n 


(2.16) 


By  Eq.  (1.34) , 


dfi  dus  dQ.  do),  .  t  dfl 

_ n  _  _ n  _ b  _  _ b  _d_X.  _ £ 

dt  dt  ’  dt  "  dt  dtlRJ  *  dt 


*-  0)  + 


(2.17) 


Introducing  Eqs.  (1.34),  (2.10),  (2.16),  and  (2.17)  into  Eq.  (2.14)  we  get 

3o)  t  t  A  2 

M„  =  Msr  •  R  “aJ  *  twb  *  P'b"  -  (il  '  l3,ua)  ‘  *1  RE 


OU).  ,  i 

Klb  "  VaT  *  5T(r)!  -  “n[i3“  -  (11  '  l3)“a1 


a  t 

Ma  =  VaT  +  R  Wn  +  w) 


(2.18) 


It  is  to  be  noted  that  the  tortuosity  1/E  enters  these  formulas  only  in  the 

final  term  of  the  equation  for  the  gyroscopic  couple  M^. 

If  the  axis  of  the  tube  is  a  plane  curve,  1/E  =  0.  The  equations  for  a 

straight  tube  are  obtained  by  setting  1/E  =  0  and  1/R  =  0.  If  the  tube  is 

immovable,  u)  =  to.  «  u>  =0.  Then  Eq.  (2.18)  reduces  to 
’  n  b  a  M 


•  •  » 

ln  =  R(l3W  11  ’  Mb  =  *1  dt^ib  ’  Ma  =  X3W 


(2.19) 
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If  the  tube  is  straight  and  rigid,  1/R  =  0  and  x  =  0.  Also,  to  is 
independent  of  s;  i.e.,  to  =  to(t) .  The  Frenet  formulas  (Eq.  (1.2))  reduce  to 

aft  6  aS  fi  at  n 

3s  "  Z  ’  3s  "  "  S  *  3s  '  U 

The  tortuosity  of  a  straight  line  is  indeterminate.  Consequently,  Z  should 
cancel  from  the  equations  for  M^,  M^,  M&.  This  condition  provides  a  partial 
check  on  the  theory. 

Equation  (1.34)  yields 

-  U)n  ,  =  (Ob  ,  -  0)a  ♦  (0  (b) 


Since  co^  =  to  •  fi,  etc.,  and  to  =  <o(t)  ,  Eq.  (a)  yields 


3to 
_ n 

3s 


3to, 
_ b 

3s 


3(0 
_ a 

3s 


0 


Consequently, 


dto 

3a) 

n 

k , 

d(o, 

b 

3u)b 

c 

£ 

dto  3u> 

a  a 

dF~  = 

at 

■f 

I  b  ’ 

dt 

3t 

•  7  w  , 
l  n 

dt  3t 

(c) 


(d) 


With  Eqs.  (b)  and  (d) ,  Eqs.  (2.10)  and  (2.14)  yield 


3<o 

Mn  =  ll  ?r  -  (il  ■  lS>“b  “a  *  V  “b 


“b- 


(i,  -  i.)to  (0  -  i.io  to 

3  n  a  3  n 


Ma  =  i3^a  *  “> 


(2.20) 


Equations  (2.20)  reduce  to  the  Euler  equations  if  to  =  0. 

It  appears  that  3a>n/3t  and  3io^/3t  depend  on  time  rates  of  change  of  n 
and  6.  However,  since  3n/3t  -  to  x  fi,  oo  •  3n/3t  =  0.  Consequently,  since 
Wn  =  to  *  n. 


n 


3w 

3t 
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Likewise, 


3tob  3w 

3t  3t 

Consequently,  3<on/3t  and  8w^/3t  do  not  depend  on  time  derivatives  of  fi  and  6. 
Therefore,  in  Eq.  (2.20),  (n,  b,  t)  may  be  any  right-handed  orthogonal  triad 

/s 

of  unit  vectors,  such  that  t  coincides  with  the  axis  of  the  tube.  Since  all 
cross  sections  of  the  tube  have  the  same  angular  velocity,  oo  is  simply  the 
angular  velocity  of  the  tube  relative  to  a  Galilean  reference  frame. 
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SECTION  3 


RESPONSE  OF  A  TAPERED  ELASTIC  CANTILEVER  GUN  TUBE  TO  EXCITATION  BY  THE 
PROJECTILE  AND  PRESCRIBED  MOTION  AT  THE  BREECH 


3.1  INTRODUCTION 

In  this  section  the  theory  in  Section  1  is  used  to  determine  the  motion 
of  a  tapered  cantilever  tube  that  is  actuated  by  the  projectile  and 
prescribed  motion  at  the  base  of  the  tube.  The  section  properties  of  the 
tube  are  arbitrary  functions  of  the  axial  coordinate  s.  Initial  droop  due 
to  gravity  is  admitted.  The  deflections  and  twist  of  the  tube  are 
represented  as  series  of  flexural  and  torsional  eigenfunctions  of  a  uniform 
cantilever  beam.  The  coefficients  in  these  series  are  functions  of  time. 
Such  series  have  the  capacity  to  converge,  in  the  least-square  sense,  to 
the  exact  solution  of  the  problem,  since  the  eigenfunctions  are  complete 
sets  of  functions.  The  coefficients  in  the  series  are  generalized 
coordinates  of  the  tube.  By  means  of  Lagrange's  equations,  they  are 
represented  as  the  solution  of  certain  coupled  non -homogeneous  ordinary 
linear  differential  equations  of  second  order  with  time-dependent  coeffi¬ 
cients. 


3.2  THE  LAGRANGIAN  FUNCTION 

The  tube  is  considered  to  be  horizontal,  and  the  y-axis  is  directed 
downward.  The  z-axis  coincides  with  the  undeflected  axis  of  the  tube. 
Accordingly,  the  potential  energy  of  the  projectile  is 

Up  -m  g  y(5,t)  (3.1) 

A  cross  section  of  the  tube  is  required  to  have  the  same  moment  of 
inertia  I  about  all  diametral  axes.  Consequently,  the  polar  moment  of 
inertia  of  a  cross  section  of  the  tube  is  21.  Accordingly,  the  kinetic 
energy  of  the  tube  is 


Ftube  ~  2  p 


s(*t2  +  yt2)ds  +  p 


I«  ds 
a 


(3.2) 


Equation  (3.2)  includes  the  torsional  kinetic  energ) ,  but  the  rotary  kinetic 
energy  due  to  the  deflections  (x,y)  has  been  neglected. 
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With  Eqs.  (1.44)  and  (1.78),  Eq.  (3.2)  yields 


tube  2 


1  2  2  &  fs 

2  p]o  3(xt  *  yt  )4s  ♦  p|o  I[*t  *  jo(xst  y 


-  v  .  x  )dsl  ds 
ss  •  st  ss'  1 


Discarding  cubic  and  quartic  terms  in  (x,  y,  ij>)  ,  we  get 


tube  2 


\  p[  S(x  2  +  y  2)ds  +  pf  I\pt2  ds 

1  Jo  x  t  j0 


(3.3) 


Since  the  twist  of  the  tube  is  t  =  the  potential  energy  of  the  tube 


is 


rZ  (Z 

utube  ■  i  J0  EI<*SS2 1  '  L  GI  ds  - <*  I  s  y 


f 

J0 


(3.4) 


The  three  expressions  in  Eq.  (3.4)  respectively  represent  the  strain  energy 
of  bending,  the  strain  energy  of  torsion,  and  the  potential  energy  due  to 
gravity. 

Equations  (1.82),  (3.1),  (3.3),  and  (3.4)  yield  the  Lagrangian  function 


L  -  T  -  U  =  i  m[xt2  .  yt2  .  2C(*s  xt  .  ys  y^) ] 


2  *l^xst  +  ^st  +  2«XS.  Xst  +  ^ss  ^st^  +  ^  ^xss  +  ^ss  ^ 


S=C 


1  . 


+  2  i3(w  +  *t  +  & sy 


s=C 


+  V 


(x  .  y  -  y  .  x  )ds 
st  ss  'st  ss' 


+  2  P 


xV  ^  a  "  « 

S(x  +  yt  ) ds  +  p  lip  ds  +  in  g  y(S,t) 
0  >0  z 


1 

-  J  E 

Jo 


l 


EI(x„2  +  yoc2)ds  -  |  Chp"2  ds  +  pg  ]  S  y  ds 


ss  '  ss 


'0 


«r 


(3.3) 


•  •  2 

The  term  (£  +  f)  has  been  omitted,  since  it  contributes  nothing  to  the 


Lagrange  equations. 


3.3  NATURAL  MODES  OF  A  CANTILEVER  BEAM 

The  n'th  natural  bending  inode  of  a  unifoi'm  cantilever  beam  (Ref.  7)  is 

f  (s)  =  cosh  3s-  cos  3  s  -  a  (sinh  3s-  sin  3  s)  (3.6) 

nK  J  n  n  n  n  ’ 

in  which  Bn&  is  the  n’th  positive  root  of  the  equation, 

cos  $£  cosh  32-  =  -1  (3.7) 

The  dimensionless  constant  a  is  defined  by 

n 

cos  3  l  +  cosh  3  £ 

°n  =  sin  VV  +"  sinh  (TT 
n  n 

Values  of  B  S,  and  a  are  given  in  Table  1.  If  n  >  5,  a  «  1  and 
n  n  b  ’  n 

3n£  ~  (2n  -  l)ir/2,  with  accuracy  at  least  to  seven  significant  figures. 


TABLE  1 

Eigenvalues  for  a  Cantilever  Beam 


n 

V 

an 

l 

1.8751041 

0.7340955 

2 

4.6940911 

1.0184664 

3 

7.8547574 

0.9992245 

4 

10.9955407 

1.0000336 

5 

14.1371684 

0.9999986 

A  few  pertinent  integrals  of  the  functions  f^(s)  are  given  below 
(Ref.  8) : 


7  D.  Young  and  R.  Felgar,  Tables  of  Characteristic  Functions  Representing 
Normal  Modes  of  Vibration  of  a  Beam ,  Engineering  Research  Series  No.  44, 
Bureau  of  Engineering  Research,  The  University  of  Texas,  Austin,  Texas, 
1949. 

^R.  P .  Felgar ,  Formulas  for  Integrals  Containing  Characteristic  Functions 
of  a  Vibrating  Beam,  Bureau  of  Engineering  Research,  Circular  No.  14,  The 
University  of  Texas,  Austin,  Texas,  1950. 
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£«(S)  f«(s)ds  = 


(  0  ,  m  j-  n 


m'  n 


(Z  /0  ,  in  ji  n 

f  "(s)  £  "(s'ds  =  ]  . 

„  m  n  1  ,,  4rt 

•'0  'B  A  ,  m  =  n 


Cf"ts)ds"2i'  l/»'(s)  v^’tv2® 


(3.9) 


Also,  derivatives  of  the  functions  f  (s)  arise.  They  are 

f  ' (s)  =  (3  [sinh  0  s  +  sin  0  s  -  a  (cosh  0s-  cos  0  s) 
n  v  ’  nl  n  n  n  n  n  ’ 


fn"($)  -  0n  [cosh  0ns  +  cos  0ns  -  an(sinh  0ns  +  sin  0n$)] 

fn'"(s)  =  0n3[sinh  0ns  -  sin  0ns  -  an(cosh  0n$  +  cos  3ns) ]  (3.10) 

The  n'th  torsional  mode  of  a  uniform  straight  tube  that  is  fixed  at  one 
end  and  free  at  the  other  is 


^n(s)  =  sin(2n  -  1)  ^ 


(3.11) 


The  following  integrals  of  these  functions  arise: 


0  ,  m  i  n 


tp  (s)i|j  (s)ds 
m  n 


J°  • 

lH/2  , 


(Z  (0 

tj)  ' (s)iii  '  (s)cls  =  j 

Jo  1,1  n  ir 


,  in  r1  n 


(2n  -  1)  ,  in  =  n 


(3.12) 


The  function  <ji(s,t)  represents  the  angular  displacement  of  a  ci'oss  section 
of  the  tube  in  its  plane.  Consequently,  the  boundary  conditions  are  iJj ( 0 , t )  = 
0,  i|i  (A,t)  =  0.  These  boundary  conditions  arc  satisfied  automatically  by 
expansion  of  <Ks,t)  in  a  truncated  series  of  the  functions  ^(s)  .  Also,  the 
boundary  conditions, 
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x(0,t)  »  x  (0, t)  =  o  ,  y(0,t)  =  yJO.t)  -  0 

S  J 

*ssC«.t)  =  XSS5(<l't)  "  °  1  W1-1’  '  W*'0  '  “ 

arc  satisfied  automatically  by  expansions  of  x(s.t)  and  y(s,t)  in  truncated 
scries  of  the  functions  f  (s) • 

3.4  EXPANSION  OP  THE  LAGRANGIAN  FUNCTION 

Prescribed  motion  of  the  breech  imposes  time -dependent  constraints  on 
the  tube.  The  Lagrange  equations  remain  valid  for  such  systems,  provided 
that  the  kinetic  and  potential  energies  are  computed  with  respect  to  a 
Galilean  reference  frame  (Ref.  9). 

As  in  Art.  1.14,  rectangular  coordinates  (x,  y,  z)  are  set  up  so  that 
the  z-axis  coincides  with  the  undeflected  and  undisplaced  axis  of  the  tube. 
The  y-axis  is  directed  downward.  Axes  (x,  y,  z)  are  attached  to  a  Galilean 
reference  frame.  The  Lagrangian  function  is  given  by  Eq.  (3.5). 

The  rectilinear  and  angular  displacements  of  the  axis  of  the  tubo  aro 
represented  as  follows: 

x  =  u(t)  +  s  <|>(t)  +  EXn(t)  fn(s) 


y  =  v (t)  -  s  0(t)  +  EYn(t)  fn(s) 


<j»  =  CCt)  +  EZn(t)  i|'n(s) 


(3.13) 


The  functions  u(t)  and  v(t)  are  the  x  and  y  components  of  displacement  at 
the  base  of  the  tube  where  s  =  0.  In  accordance  with  Eq.  (1.71),  the 
z-component  of  displacement  at  the  base  of  the  tube  is  f(t),  but  this  term 
is  irrelevant..  The  x  and  y  components  of  rotation  at  the  base  of  the  tube 
are  6(t)  and  <|>(t) ,  respectively.  The  z-component  of  rotation  at.  the  base  of 
the  tube  is  5(t) .  The  functions  f  (s)  and  *n(s)  arc  defined  by  Eqs.  (3.6), 
(3.7),  (3.8),  and  (3.11).  The  functions  u(t),  v(t.)  ,  <|>(t) ,  9(t),  and  C(t) 
are  considered  to  be  given.  The  functions  *^(t)  ,  Yn(t),  z-n(t)  are 
generalized  coordinates  of  the  tube.  The  range  of  the  subscript  a  in 


9 !!.  L.  Langhaar,  Energy  Methods  in  Applied  Mechanics,  John  Wiley  :■  Sons, 
hew  York /  1962 ,  Art.  7 -A. 
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Gq.  (3.13)  is  1,  2,  3,  ...  ,  N,  but  the  number  N  is  unspecified,  and,  for 

simplicity,  the  range  of  n  is  not  indicated  on  the  summation  signs. 

The  Lagrangian  function  is  obtained  by  substituting  Eq.  (3.13)  into 

Gq.  (3.3),  Terms  that  do  not  contain  the  dependent  variables  X  .  Y  Z  or 
.....  .  »  n '  n 

tlieir  derivatives  are  omitted  from  the  Lagrangian  function,  since  they 

cancel  from  the  Lagrange  equations.  The  result  is 

l  -  i  ,nmim  i„  *  V|u  \,CJX)  f„(E)  *  xn  *  V1#  yj yo  yr, 

*  I  liES<*.  \  *  vn)f,„'(E)  V  (E)  .  Xn 

'  Ym  VV®  V®  *  I  X„  +  Y„  W<E)  V  (E) 

*  1  K  ♦.(O  V5>  -  I  i/ErZ»  Zn  V® 

*  l3“d„  V°  *  ♦»’«>  *  1&EZ„  z„  V®  *n(E) 


i3“E£I(X.  Vn  -  Y»  V  „  V  f„"  ds)  +  ?  <>«[<**  X„ 


*»  ?»>j*  3  <„  f„  ds!  *  iS  I  *m  Us] 


S ]  +  In  g  Z  Y  f  (5) 
n 


?  ££[(Xm  X„  *  Ym  vf#  E'  V  V  dsl  '  EEtZ„  zj*  Cl  yds] 


+  PgZ[Y 


rl 


I  S  f  ds]  +  m(u  +  fi)SX  f  '  (£)  +  m(v  -  f  (£) 

1  o  11  11  >1  11  ti 

+  m£ (u  +  C*)>:xn  fn'(5)  +  m^XXn  fn(?)  -  n|0EYn  fn(0 
-  mC(v  -  cemn  fn’(C)  +  ij  «fcxn  fn'(C)  -  ij  0EYn  fn'(C) 

♦  i,WXn  fn"K)  -  i,Efe  vn  fn-.  (?)  ♦  ijiij  ♦  (O 
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+  ijttEZn  <^’(?)  ♦  i3wl(0  Xn  +  4,  Yn)fn'(C)  +  P^(u  Xn 


rJc 

S  f  ds  +■  pE  (d>  X  -  0  Y  ) 

J0  n  n  n  J 

.  . 

S  s  f  ds  +  2p£EZ 

1 0  n  nj 

I  tp  ds 
n 

(3.14) 


If  the  tube  is  uniform,  the  integrals  over  the  range  (0,&)  in  Eq.  (3.14)  can 
be  evaluated  by  means  of  Eqs.  (3.9)  and  (3.12). 

In  the  present  case,  ?  and  w  are  regarded  as  given  functions. 
Consequently  6?  and  Sy  are  zero.  Hence,  by  the  argument  in  Art.  1.10, 

<$W  =  0.  Therefore,  the  Lagrange  equations  are 


_d_  9L_  9L_ 

dt  -  ‘  axr 


d  3L 


_d_  9L_ 

dt  az 


(3.15) 


3.5  LAGRANGE'S  EQUATIONS 

Since  fn(?)  is  a  function  of  t,  Eqs.  (3.14)  and  (3.15)  yield 


E  [®f  (?)  f  (?)  +  i,  f  '  (?)  f  '(?)  +  p 
n  r  n  1  r  n 


S  f  f  ds] X 
r  n  1  n 


+  2 ?E  [m  f  (?)  f  '(?)  +  i  f  »(?)  f  "(?)]X 
n  r  ti  1  r  n  n 


+  Z  [m?  f  (?)  f  '(?)  +  m?2  f  '(?)  f  '(?)  +  m?2  f  (?)  f  "(?) 
n  r  n  r  n  r  n 


+  i.?  f  '(?)  f  "(?)  +  i.?2  f  '(?)  f  '"(?)  + 
1  r  n  v  }  1*  r  vw  n 


El  £  "  f  "  ds] X 

Jo  r  n  n 


+  i_w  f  '  (?)  E  f  '  (?)Y  +  i_Z  [w?  f  '(?)  f  "(?)  +  u) 

3  r  v  ’  n  n  3n  r  n 


f  '  f  "  ds] Y 
r  n  J  n 


=  t -mu  -  m?$  -  m?cf)  -  2  m?4>]fr(?)  +  |-m?2<j>  -  +  i^O]^’  (?) 


fX, 

-  pii  S  f  ds  -  p>j> 


S  s  f  j  s 

r 


(3.16) 


57 


■I 

E  [m  f  (0  f  ©  +  i.  f  '(O  f  '(Q  +  P  S  f  f 
n  1  r  n  1  r  n  A  r  n 

(1 


*  2&  [m  £r(«  fn'(5)  +  ix  fr'(€)  V(C)]Yn  +  n  fr(5)  V  (<0 

♦  mC2  fr'(0  fn*(0  +  mt2  fr(0  fn"(C)  +  ij?  f.'(0  f/’  (O 


*  ixr  frf(c)  fn,M(c)  + 


(S. 

El  f  »  f  "  ds] Y  -  i_w  f  '(£)!  £  •  (OX 
Jn  r  n  J  n  3  r  v  'n  n  v  ”  n 


-  i  l  to£  f  •  (CD  f  "(O  +  w 

3n  r  n 


.r« 

u  f  '  f  "  ds]X 

J0  r  n  n 


[-mv  +  2mC0  +  m£0  +  m£0  +  m  g]  £  (?)  +  [m5  0  +  +  i  u«t>l  fr •  (5) 


pv  S  f  ds  +  p0  S  s  f  ds  +  pg  S  f  ds 
Jn  r  n  r  Jn  r 


(3.17) 


„  [i3  V«>  +  H0  1  *r  *n  ds^n  +  2^3  V(«*n 


+  2  [i3e  <l>r(£)  V(5)  +  *3^  V(Q  +  2  GI  V  V  ds]Zn 


r* 

2  GI  U>  *  i/>  * 
J0  vr  •! 


■i3(w  +  0^r(Q  -  2pcJ  I  t|>r  ds 


(3.18) 


Equations  (3.16),  (3.17),  and  (3.18)  are  ordinary  linear  non- 
homogeneous  differential  equations  of  second  order.  They  have  the  following 
form: 

Z  X  F  !(t)  +  Z  X  F  2(t)  +  Z  X  F  3(t)  +2  Y  F  4(t) 

n  n  rn  v  J  n  n  rn  v  n  n  rn  n  n  rn 

+  Z  Y  F  5(t)  =  K  X(t) 
n  n  rn  r 
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Z  Y  F  X(t)  +  E  Y  F  2(t)  +  E  Y  F  3(t)  -  E  X  F  4(t) 
n  n  rn  '  7  n  n  rn  v  n  n  rn  v  1  n  n  rn  ’ 


E  X  F  5(t)  K  2(t) 
n  n  rn  ^  ;  r  v  J 


E 

n 


Z  H  „A(t)  +  E  Z  H  2(t)  +  E  Z  H  3(t)  -  K  3(t) 
n  rn  K  n  n  rn  v  J  n  n  rn  v  '  r  v  ' 


(3.19) 


Since  £  and  in  are  regarded  as  known  functions  of  t,  the  coefficients  F  \ 

H  **  >  K  are  presumably  known,  at  least,  in  tabular  form.  Although  they 

are  rather  complicated,  they  can  be  programmed  for  a  computer.  It  is 

noteworthy  that  the  root  excitation  functions  u(t) ,  v(t) ,  9(t),  <}>(t),  C(t) 

do  not  enter  into  the  functions  F  ^  or  H  •* :  they  affect  only  the 

rn  rn  ’  7  ' 

functions  K  ■* .  The  case  of  a  canti  lever  tube  that  is  fixed  at  the  root  is 

r 

obtained  by  setting  u  =  v  =  0  and  0  =  <j>  =  5  *  0. 

If  the  gun  is  initially  at  rest,  the  initial  conditions  are 


X  (0)  =  Y  (0)  =  Z  (0)  =  0  and  X  (0)  =  Y  (0)  =  Z  (0)  =  0  (3.20) 

n  n  n  n  n  n  ' 


Equations  (3.19)  and  (3.2G)  present  an  initial-value  problem  of  a  type  for 
which  numerical  methods  are  available.  The  fact  that  the  Z r  equations  are 
separated  from  the  others  is  helpful.  However,  there  is  coupling  between 
the  X^  and  Y^  equations.  Accordingly,  a  vertical  oscillation  of  the  tube 
excites  a  horizontal  oscillation.  This  coupling  vanishes  if  the  spin  w  is 
zero.  The  coupling  manifests  gyroscopic  action  of  the  projectile. 
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SECTION  4 


GYROSCOPIC  ACTION  OF  A  BALANCED  SPINNING  PROJECTILE 
IN  A  MOVABLE  RIGID  CURVED  TUBE 


4.1  INTRODUCTION 

The  barrel  of  a  gun  is  unavoidably  slightly  curved  because  of  effects 
of  gravity,  temperature  gradients,  manufacturing  imperfections,  etc.  The 
spinning  projectile  consequently  exerts  a  gyroscopic  couple  that  tends  to 
bend  the  tube  sideways.  As  Section  3  shows,  the  analysis  of  this  action  for 
an  actual  gun  is  complicated,  although  it  appears  to  lie  within  the  scope  of 
numerical  methods  for  differential  equations.  However,  insight  is  gained 
by  studying  much  simpler  problems  that  are  not  without  practical  significance. 
Consequently,  in  this  section,  the  motion  of  a  rigid  tube  whose  axis  is  a 
plane  curve  is  analyzed  for  a  gun  that  is  hinged  at  the  breech  so  that  the 
tube  can  swing  sideways.  A  resisting  moment  M  that  depends  arbitrarily  on 
the  side  sway  <J)  and  its  time  derivative  <{>  is  introduced.  The  moment  M  that 
yields  <f  =  0  is  that  which  is  experienced  by  a  rigid  immovable  gun. 

The  motion  of  the  projectile  in  the  tube  is  considered  to  be  prescribed. 
Balloting  is  disregarded.  The  projectile  is  considered  to  be  perfectly 
balanced.  The  geometric  axis  of  the  projectile  is  accordingly  tangent  to 
the  axis  of  the  tube. 

For  comparative  purposes,  two  different  methods  of  solution  are 
employed.  The  first  treatment  is  based  on  the  principle  of  angular  momentum. 
The  second  treatment  is  based  on  Lagrange's  equation. 

4.2  LATERAL  MOTION  OF  A  HINGED,  RIGID,  CURVED  TUBE 

The  vertical  axis  a-a  (Figure  6)  is  taken  to  be  a  hinge  line.  The 
hinge  contains  a  spring  and  a  damper,  which  may  be  nonlinear.  The  hinge 
allows  the  tube  to  swing  horizontally.  The  angular  displacement  of  the  tube 
about  the  hinge  is  denoted  by  <j>.  The  angular  velocity  of  the  tube  is  <)>,, 
where  the  dot  denotes  the  time  derivative.  The  axis  of  the  tube  is 
considered  to  be  a  plane  curve  with  curvature  1/R.  The  spin  of  the 
projectile  relative  to  the  tube  is  w(t) . 

The  angular  velocity  of  the  projectile  has  components  <p  and  to  in  the 

plane  of  the  axis  of  the  ’•ube.  The  axial  and  normal  components  of  vector 
•  •  •  * 

<J>  are  -<f>  sin  0  and  <p  cos  0,  as  shown  by  Figure  6.  The  axial  component  of  <p 
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detracts  from  the  spin  vector,  so  the  net  absolute  axial  component  of 
angular  velocity  of  the  projectile  is  w  -  <J>  sin  0.  The  transverse  component 
of  angular  velocity  of  the  projectile  is  £/R,  where  £  is  the  speed  of  the 
projectile  relative  to  the  tube.  Accordingly,  the  absolute  angular  velocity 
components  of  the  projectile  are 

•  f*  • 

-ft  =  cos  0,  ft^  =  =  (*)-<{>  sin  0  (4.1) 

Consequently,  the  principal  components  of  angular  momentum  of  the  projectile 
about  axes  through  the  center  of  mass  of  the  projectile  are 

» 

Hj  ~  ij  4-  cos  6  ,  1^2  =  *1  R  ’  =  i3^°  "  sin  (4.2) 

where  i^  is  the  moment  of  inertia  of  the  projectile  about  the  transverse 
axis  through  its  center  of  mass,  and  i^  is  the  moment  of  inertia  of  the 
projectile  about  its  longitudinal  axis.  The  angular  momentum  of  the 
projectile  about  the  axis  of  the  hinge  is 

Hj  cos  0  -  H3  sin  0  +  m  x2  <j>  (4.3) 

where  m  is  the  mass  of  the  projectile  and  x  is  the  distance  from  the  center 
of  mass  of  the  projectile  to  the  axis  of  the  hinge.  Equation  (4.3)  must  be 
augmented  by  the  angular  momentum  J(|>  of  the  tube  and  the  breech,  where  J  is 
the  moment  of  inertia  of  the  tube  and  the  attached  rotating  part  of  the 
breech  about  the  axis  of  the  hinge,  Consequently,  if  the  angular  momentum 
of  the  charge  is  disregarded,*  the  angular  momentum  of  the  system  about  the 
hinge  line  is 

H  =  J<|>  +  Hj  cos  9  -  Hj  sin  9  m  (4.4) 

Consequently,  by  Eq.  (4.2), 

2  2  2  • 

H  =  (J  +  ij  cos  0  +  i^  sin  0  +  in  x  )<j>  -  sin  9  (4.S) 

;t'  ' 

The  angular  momentum  of  the  charge  may  be  introduced  in  an  empirical  way 
by  augmenting  the  mass  of  the  projectile  by  a  part  of  the  mass  of  the 


The  angular-momentum  principle  is  expressed  by  the  equation 


(4.6) 


where  M(<fi  ,<£)  is  the  resisting  moment  of  the  spring  and  damper  in  the  hinge. 
The  quantity  dH/dt  is  the  substantial  derivative;  i.e.,  it  is  the  time  rate 
of  change  of  H  with  due  regard  for  the  time  dependence  of  0,  x,  and  w.  It 
can  be  seen  that  d0/dt  =  £/Rt  doi/dt  =  u>,  and  dx/dt.  =  £  cos  0.  Consequently, 
Eqs.  (4.5)  and  (4.6)  yield 

2  2  2  “  £ 

(J  +  ij  cos  0  +■  i^  sin  0  +  m  x  )<j>  +  2[(ig  -  ij)  ^  sin  0  cos  0 


+  in  £  x  cos  0]4>  -  igio  sin  0  -  i^  j*-  uj  cos  0  +  M(<j>,<f>)  =  0 


(4.7) 


Equation  (4.7)  is  a  second-order  ordinary  differential  equation  that 
determines  tj>(t)  ,  if  the  initial  values  <f>(0)  and  tj>(0)  are  given.  The 
function  M((f>,<t0  may  be  nonlinear.  Otherwise,  Eq.  (4.7)  is  linear. 

If  M  =  0,  the  gun  swings  freely.  Then,  Ea.  (4.6)  yields  H  =  constant. 
Therefore,  if  M  =  0, 


2  2  2  ♦ 

(J  +  ij  cos  0  +  i3  sin  0  +  m  x  )<f>  -  i^u)  sin  0 
2  2  2  » 

=  (J  +  ix  cos  eQ  +  i3  sin  0Q  +  m  xQ  )<|>0  -  sin  0Q  (4.8) 


•  ♦ 

where  0Q,  x^,  are  initial  values.  Equation  (4.8)  determines  (j> 

explicitly. 

Since  0  ordinarily  is  a  very  small  angle,  it  is  reasonable  to  make  the 
approximations  sin  9  ~  0  and  cos  0  ~  1  in  Eqs,  (4,7)  and  (4.8).  If  these 
approximations  are  made,  and  if  <f>Q  =  8q  =  0,  Eq.  (4.8)  yields 


<Ht) 


rt 


w  0  dt 


0  J  +  ij  +  m  x 


(4.9) 
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4.3  LATERAL  MOTION  DERIVED  FROM  LAGRANGE’S  EQUATION 

The  velocity  components  of  the  center  of  mass  of  the  projectile 
(Figure  6)  are 

V  =  0  ,  v2  =  x  |  ,  v3  =  £ 


where  E,  is  the  distance  that  the  projectile  has  moved  along  the  axis  of  the 
tube.  Consequently ,  the  kinetic  energy  of  translation  of  the  projectile  is 


C2) 


The  principal  rotation  components  of  the  projectile  are  given  by  Eq.  (4.1j. 
Consequently,  the  kinetic  energy  of  rotation  of  the  projectile  is 

•  2 

1  *2  2  E  1  •  2 

y  ij(<^  cos  0  +  —■)  +  j  i3(o)  -  |  sin  0j 

R 

Accordingly,  the  kinetic  energy  of  the  system  is 

*  2 

1  *2  1  2*2  9 2  1  *2  2  E 

T  =  j  J  <r  +  ~2  ">(x  4>  +  C)  +  |  ixC4>  cosZ  0  +  *j) 

R 

+  |  i3(w  -  <f>  sin  0) 2  (4.10) 

The  virtual  work  of  the  external  forces  is 


<5W  =  -M<$<j>  =  Q<S<p 


Accordingly,  the  generalised  external  force  is  Q  =  -M.  Lagrange's 
equation  for  <J>  is 


_d_  3T  _  3T  n 

dt  9r 9* 


(4.11) 


In  the  present  case,  3T/94>  =0.  By  Eq.  (4.10), 


ri’T  2  2  ?  • 

—  =  (J  +  i,  cos  0  +  i  sin  9  +  m  x  )<{>  -  i _oo  sin  0  (4.12) 

34)  1  3  3 
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Consequently,  in  view  of  Eq.  (4.5),  3T/B$  is  identified  as  the  angular 
momentum  H  of  the  system  about  the  hinge.  Therefore,  Eq.  (4.11)  is 
identical  to  Eq.  (4.6),  and  Lagrange’s  equation  leads  to  Eq.  (4.7). 


4.4  MOMENTS  IN  A  RIGID  IMMOVABLE  CURVED  TUBE 
If  <j>  =  0,  Eq.  (4.7)  gives 

M  =  i3(w  sin  0  +  ~  cos  8)  (4.13) 

Equation  (4.13)  gives  the  sidewise  moment  on  the  tube  at  the  breech,  if  the 
tube  is  immovable.  It  acts  to  oppose 
The  rifling  torque  is  (Figure  7) 

i  w  =  M  =  M 

j  r 


This  is  the  torque  exerted  on  the  projectile  by  the  tube.  The  torque 
exerted  on  the  tube  by  the  projectile  is  -M_.  The  component  of  this  torque 
on  the  axis  of  the  hinge  is 

i^uj  sin  0 


The  driving  force  of  the  gases  exerts  no  moment  about  the  hinge.  There  is 
no  force  on  the  projectile  transverse  to  the  plane  of  the  axis  of  the  tube. 
Consequently,  the  reaction  of  the  forces  on  the  projectile  exerts  no 
moment  about  the  hinge  line.  The  gyroscopic  couple  that  the  projectile 
exerts  on  the  tube  (Figure  7)  is  denoted  by  M^.  Equilibrium  of  moments 
about  the  hinge  line  yields 

-M  +  M  cos  0  +  M_  sin  6=0 
g  3 


Therefore,  by  Eq.  (4.13), 


-i3(u>  sin  0  + 


cos  0)  +  M  cos  0  +  i,w  sin  0=0 
g  3 


This  reduces  to 


M 


g 


(4.14) 
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■.  "w  A.  ^  A’  *. 


M3H,W  =  Mr 


Figure  7.  Components  of  Moments  on  a  Projectile  in  a  Fixed  Rigid  Curved 
Tube 


The  positive  sense  of  M  is  indicated  by  Figure  7,  if  the  spin  of  the 

8 

projectile  is  that  of  a  right-hand  screw  advancing  along  the  tube.  Although 
1/R  is  very  small,  w  and  v  are  very  large.  Consequently,  Eq.  (4.14) 
indicates  that  the  gyroscopic  action  of  the  projectile  might  bend  the  tube 
appreciably. 

The  moment  component  about  the  binormal  is  =  ij  0.  Hence, 


=  i  —(£)  =  i  [!.  +  £2  (i)l 
1  dtV  llR  *  dslR,j 


(4. IF) 


where  s  is  arc  length  on  the  axis  of  the  tube. 

If  the  axis  of  the  rigid  tube  is  a  space  curve,  it  can  be  shown  that 
Eq.  (4.14)  is  generalized  as  follows: 


M  =  1  — —  -  ]  -2 — 

g  3  R  1  RE 


(4.16) 
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4.5  MOVING  RIGID  TUBE  AS  A  SPECIAL  CASE  OF  A  FLEXIBLE  TUBE 

As  an  example,  Eq.  (1.41)  is  applied  to  the  problem  treated  in  this 
section.  The  axis  of  the  tube  is  a  plane  curve.  The  tube  is  rigid,  and 
the  vertical  axis  (Figure  8)  is  a  hinge  line.  As  before,  the  angular 
displacement  of  the  tube  about  the  hinge  is  *#» • 


0 


Figure  8.  Rotating  Rigid  Tube 

Since  the  tube  is  rigid,  the  rate  of  twist  is  zero;  i.e.,  9x/3t.  =  0. 
The  binormal  S  is  perpendicular  to  the  plane  of  Figure  8.  It  is  directed 
away  from  the  reader.  The  vector  3t/8t  has  the  direction  of  b.  Its 
magnitude  is  4>  cos  0,  Accordingly,  S  *  3t/9t  =  <j>  cos  0,  and  Eq.  (1.41) 
yields 
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*  jf**»  i.  >  * 


I 


* 

JT  *  it  *  cos  8  -  0 


Also,  1/R  =  30/ 9s.  Consequently, 


a  I  oc  ...  a 
__  =  _  cos  8 


Integration  yields 


toa  =  -(ji  sin  0  +  f(t) 


(4.17) 


It  has  been  shown  in  Art..  4.2  that  f(t)  -  0,  but  there  seems  to  be  nothing 
in  the  theory  of  the  flexible  tube  that  determines  f(t).  With  f(t)  =  0, 
Eq.  (1.33)  and  Eq.  (4.17)  yield  the  angular  velocity  components  of  the 
pi'ojectile : 

• 

♦  r  ♦ 

Q,  =  -<i>  cos  0,  ft.  =  £  ,  Q  =  u>  -  <{>  sin  0 

H  UK  a 

These  results  agree  with  Eq.  (4.1). 
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SUCTION  5 
CONCLUSIONS 


Vector  and  scalar  formulas  for  the  velocity,  the  acceleration,  tho 
angular  velocity,  and  the  kinetic  energy  of  a  geometrically  perfect 
projectile  in  a  concentric  flexible  tube  are  derived  rigorously  in  Section  1. 
Approximations  of  these  formulas  for  an  initially  straight  tube  also  are 
developed.  The  relation  between  twist  of  tho  tube  and  angular  velocities 
of  cross  sections  of  the  tube  is  complicated  by  curvature  of  the  axis  of 
the  tube.  This  matter  is  examined  in  Art-.  1.9.  Also,  evaluation  of  the 
virtual  work  of  the  forces  associated  with  the  projectile  is  primarily  a 
kinematic  problem.  It  is  investigated  in  Art.  1.10.  Aside  from  the  complex 
phenomenon  of  balloting,  Section  l  lays  a  rigorous  kinematic  foundation  for 
gun  dynamics.  Section  2  deals  with  the  forces  and  moments  acting  on  a 
dynamically  balanced  projectile  in  a  flexible  tube.  It  provides  formulas 
for  the  rifling  torque  and  the  gyroscopic  couple. 

Section  3  treats  the  action  of  a  spinning  pi'ojectile  on  a  tapered 

elastic,  cantilever  tube  that  has  prescribed  motion  at  the  breech.  The 

deflections  and  twist  of  the  tube  are  expanded  in  series  of  natural  bonding 

modes  and  torsional  modes  of  a  uniform  cantilever  beam.  Since  these  modes 

constitute  complete  sets  of  functions,  truncated  series  of  them  can 

represent  tho  deflections  and  twist  of  the  tube  to  any  desired  degree  of 

accuracy  (in  tho  least-square  sense),  irrespective  of  variable  taper  of  the 

tube.  The  coefficients  (X  ,  Y  ,  Z  )  in  tho  modal  series  are  time-dependent 

n  n  n 

generalized  coordinates  for  the  tube.  Lagrange's  equations  provide  linear, 
second-order,  non-homogeneous,  oi-dinaty,  differential  equations  for  (X^,  Y^, 
Zn)  .  Although  the  coefficients  in  these  differential  equations  are 
complicated  functions  of  time,  the  differential  equations  may  be  expected 
r.o  be  amenable  to  numerical  methods  that  can  bo  programmed  for  a  digital 
computer.  Although  the  theory  in  Section  3  is  not  immediately  applicable  to 
a  gun  for  which  the  motion  of  the  breech  is  unknown,  it  may  be  assumed 
tentatively  that  the  motion  of  the  breech  complies  approximately  with  that 
of  a  completely  rigid  gun  with  the  type  of  support  and  recoil  mechanism 
that  is  under  consideration  (Ref.  1). 

Although  gyroscopic  action  of  the  projectile  is  inherent  in  the 
behavior  of  a  flexible  tube,  the  phenomenon  is  combined  with  centrifugal 
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effects  of  the  projectile  in  the  deflected  tube  and  other  complications. 
Section  4  isolates  the  phenomenon  in  a  setting  that  has  practical  elements. 
A  rigid  curved  tube  whose  axis  lies  in  a  vertical  plane  is  hinged  at  the 
breech,  so  that  the  tube  can  swing  sideways.  A  resisting  moment  M,  that 

* 

deponds  arbitrarily  on  the  angle  <f>  of  side  sway  and  its  time  derivative  <j>, 
is  introduced.  The  moment  M  that  yields  <j>  =  0  is  that  which  is  experienced 
by  a  rigid  immovable  gun.  A  single  second-order  ordinary  differential 
equation  that  determines  4> (t)  is  derived.  The  function  may  be 

nonlinear.  Otherwise,  the  differential  equation  is  linear.  A  numerical 
study  of  the  differential  equation  should  be  instructive. 

Gyroscopic  couples  acting  on  the  projectile  are  reacted  on  the  tube. 
According  to  elemontary  beam  theory,  the  curvature  of  the  tube  is 
proportional  to  the  bending  moment.  Consequently,  a  gyroscopic  couple 
(ideally  conceived  to  act  at  the  cross  section  of  the  tube  where  the  center 
of  mass  of  the  projectile  lies)  introduces  a  stepwise  discontinuity  in  the 
curvature.  This  anomaly  portends  a  puzzling  mathematical  question,  since 
the  gyroscopic  couple  depends  on  the  local  curvature  of  the  tube,  but  the 
curvature  is  indeterminate  at  the  point  where  the  couple  is  conceived  to 
act.  However,  numerical  methods  tend  to  smooth  over  discontinuities.  For 
example,  any  linear  combination  of  natural  inodes  of  the  tube  is  continuous. 
Likewise,  a  piecewise  polynomial  is  continuous  with  all  its  derivatives, 
except  at  the  junctions  of  the  polynomial  segments. 
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US  Army  Tank  Automotive  Research 
and  Development  Command 
ATTN:  DRDTA-RH,  Dr.  IV. F.  Banks 
DRDTA,  Dr.  E.  Patrick 
Dr.  Jack  Parks 
Warren,  MI  48090 

1  Director 

US  Army  TRAD0C  Systems 
Analysis  Activity 
ATTN:  ATAA-SL,  Tech  Lib 
White  Sands  Missile  Range,  NM  88002 

2  President 

US  Army  Armor  and  Engineer  Board 
ATTN:  ATZK-AE-CV 

ATZK-AE-IN,  Mr.  Larry  Smith 
Fort  Knox,  KY  40121 

2  Commander 

US  Army  Research  Office 
ATTN:  COL  L.  Mittenthal 
Dr.  E.  Saibel 
P.O.  Box  12211 
Research  Triangle  Park 
NC  27709 

3  Commander 

US  Army  Research  Office 

P.O.  Box  12211 

ATTN:  Technical  Director 

Engineering  Division 
Metallurgy  fi  Materials 
Division 

Research  Triangle  Park,  NC  27709 

1  Commander 

US  Army  Research  Office 

ATTN:  Dr.  J.  Chandra 

Research  Triangle  Park,  NC  27709 
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2  Pro  j  ect  Manager 
Nuclear  Munitions 
ATTN :  DRCPM-NUC 
Dover,  NJ  07801 

2  Project  Manager 

Tank  Main  Armament  Systems 
ATTN :  DRCPM-TMA 
Dover,  NJ  07801 

2  Pro j  ect  Manager 

Division  Air  Defense  Gun 
ATTN :  DRCPM-ADG 
Dover,  NJ  07801 

1  Product  Manager  for  30mm  Ammo. 
ATTN:  DRCPM-AAH~30mm 

Dover,  NJ  07801 

2  Product  Manager 

Ml  10E2  Weapon  System,  DARC0M 
ATTN:  DRCPM-M1 10E2 
Rock  Island,  IL  61299 

4  Director 

US  Army  Mechanics  and 
Materials  Research  Center 
ATTN:  Director  (3  cys) 
DRXMR-ATL  (1  cy) 
Watertown,  MA  02172 

2  Commander 

US  Army  Materials  and 
Mechanics  Research  Center 
ATTN:  J.  Mescal  1 

Tech.  Library 
Watertown,  MA  02172 

1  Commander 

US  Army  Training  and 
Doctrine  Command 
ATTN:  TRAD0C  Lib,  Mrs.  Thomas 
Fort  Monroe,  VA  23651 
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1  Commander 

US  Army  Air  Defense  Center 
Fort  Bliss,  TX  79916 


1  Commander 

Naval  Research  Laboratory 
Washington,  DC  20375 

1  Commander 

Naval  Sea  Systems  Command 
Washington,  DC  20362 

2  Commander 

Naval  Sea  Systems  Command 
ATTN:  SEA-62R,  John  W.  Murrin 
SEA-624B,  John  Carroll 
Washington,  DC  20362 

2  Commander 

Naval  Sea  Systems  Command 
(SEA- 035 13) 

ATTN :  L.  Pasiuk 
Washington,  DC  20362 

1  Commander 

Naval  Air  Systems  Command 
ATTN:  AIR-604 
Washington,  DC  20360 

1  Commander 

Naval  Ship  Engineering  Center 
Washington,  DC  20362 

1  Superintendent 

Naval  Postgraduate  School 
ATTN:  Dir  of  Lib 
Monterey,  CA  93940 

1  Commander 

Naval  Air  Development  Center 
Johnsville 

Warminster,  PA  18974 
1  Commander 

Naval  Research  Laboratory 
Washington,  DC  20375 
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1  Commander 

David  W.  Taylor  IJIaval  Ship 
Rsch  S  Development  Center 
Betliesda,  MD  20084 

3  Commander 

Naval  Research  Laboratory 
ATTN:  Mr.  IV. J.  Ferguson 
Dr.  C.  Sanday 
Dr.  H.  Pusey 
Washington,  DC  20375 


1  Commander 

Naval  Surface  Weapons  Center 
ATTN:  G-13,  W.D.  Ralph 
Dahlgren,  VA  22448 

6  Commander 

Naval  Surface  Weapons  Center 
ATTN:  Code  X21,  Lib 
E,  Zimet,  R13 

R. R.  Bernecker,  R13 

J. IV.  Forbes,  R13 

S. J.  Jacobs,  RIO 

K.  Kim,  R13 

Silver  Spring,  MD  20910 


1  Commander 

Naval  Weapons  Center 
China  Lake,  CA  93555 

4  Commander 

Naval  Weapons  Center 
ATTN:  J.  O'Malley 
D.  Potts  . 

Code  3835,  R.  Sewell 
Code  3431,  Tech  Lib 
China  Lake,  CA  93555 

1  Commander 

Marine  Corps  Development 

and  Education  Command  (MCDEC) 
ATTN:  Class  Ctl  Ctr 
Quantico,  VA  22134 

3  Commander 

Naval  Weapons  Center 
ATTN:  Code  4057 

Code  401 lj  B.  Lund strom 
Code  5835  M,  Backman 
China  Lake,  CA  93555 

1  Commander 

Naval  Ordnance  Station 
Indian  Head,  MD  20640 


3  Commander 

Naval  Surface  Weapons  Center 
ATTN:  Code  E-31,  R.C.  Reed 

M.T.  W'alchak 
Code  V-14,  W.M.  Hinckley 
Silver  Spring,  MD  20910 

1  Commander 

Naval  Surface  Weapons  Center 
Silver  Spring,  MD  20910 

5  Commander 

Naval  Surface  Weapons  Center 
ATTN:  Code  G-33,  T.N.  Tschirn 
Code  N-43,  J.J.  Yagla 

L.  Anderson 
G.  Soo  Hoo 

Code  TX,  Dr.  IV. G.  Soper 
Dahlgren,  VA  22448 


2  Commander 

Naval  Ordnance  Station 
ATTN:  Code  5034,  Ch,  Irish,  Jr 

T.C.  Smith 

Indian  Head,  MD  20640 

1  Office  of  Naval  Research 

ATTN:  Code  ONR  439,  N.  Perrone 
Department  of  the  Navy 
800  North  Quincy  Street 
Arlington,  VA  22217 
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2  AFRPL 

ATTN:  IV.  Andrepont 
T.  Park  . 

Edwards  AFB,  CA  95523 
1  AFOSR 

Bolling  AFB,  DC  20332 


3  AFATL  (DLA) 

ATTN:  K.  Dittrich;  DLJM 
D.  Davis;  DLDL 
Eglin  AFB,  FL  32542 

1  ADTC/DLODL,  Tech  Lib 
Eglin  AFB,  FL  52542 

1  AFWL/SUL 

Kirt land  AFB,  NM  8711S 


1  New  Mexico  Institute  of 
Mining  and  Technology 
Terra  Group 
Socorro,  NM  87801 

1  Director  3 

Los  Alamos  Scientific  Lab 
P.0.  Box  1663 
Los  Alamos,  NM  87544 

1  AFML  (LLN/Dr.  T.  Nicholas) 
bright -Patterson  AFB,  OH  45433 

2  ASD  (XROT,  Gerald  Bennett) 

ENFTV,  Martin  Lentz 
Wright-Patterson  AFB, 

OH  45433 


1  U.S.  Department  of  the  Interior  ^ 

Bureau  of  Mines  i 

Pittsburgh  Technical  Support  j 

Center  • 

ATTN:  Dr.  S.G.  Sawyer  J 

Pittsburgh,  PA  15213  • 

2  Battel le  Pacific  Northwest  Laboratories 

ATTN:  Dr.  F.  Simonen  j 

Mr.  E.M.  Patton  • 

P.O.  Box  999  i 

Richland,  KA  99352  j 

i 

1  Director  j 

Lawrence  Livei'more  Laboratory  5 

P.O.  Box  808  < 

Livermore,  CA  94550  j 

1  Director  ) 

Lawrence  Livermore  Laboratory  j 

ATTN:  D.  Burton,  L200  3 

P.O.  Box  808  ] 

Livermore,  CA  94550  j 

1  Director  j 

Lawrence  Livermore  Laboratory  j 

ATTN:  J.  Fleck,  L71  (Mail  Code) 

P.O.  Box  808  j 

Livermore,  CA  94550  | 

i 

1  Director  ] 

Lawrence  Livermore  Laboratory  | 

ATTN:  E.  Farley,  L9  (Mail  Code)  j 

P.O.  Box  80S  3 

Livermore,  CA  94550  j 


3  Director 

Lawrence  Livermore  Laboratory 
ATTN:  Dr.  R.H.  Toland,  L-424 
Dr,  M.L.  Wilkins 
Dr .  R .  Werne 
Livermore,  CA  94550 
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2  Aerospace  Corporation 
ATTN':  Mr.  L.  Rubin 
Mr.  L.  G.  King 
2350  E.  El  Segundo  Boulevard 
El  Segundo,  CA  90245 

1  Headquarters 

National  Aeronautics  and 
Space  Administration 
Washington,  DC  20546 

1  Sandia  Laboratories 
ATTN:  M.L.  Merritt 
P.0.  Box  5800 
Albuquerque,  NM  87115 

1  Director 

Jet  Propulsion  Laboratory 
ATTN :  Lib  (TD) 

4800  Oak  Grove  Drive 
Pasadena,  CA  91103 

1  H.P.  White  Laboratory 
Bel  Air,  MD  21014 

1  DuPont  Experimental  Labs 
Wilmington, 

DE  19801 

1  Materials  Research  Laboratory,  Inc 
1  Science  Road 
Glenwood,  IL  60427 

1  Princeton  Combustion  Research 

Laboratories,  Inc. 

ATTN:  Prof.  M.  Summerfield,  Pres. 
1041  U.S.  Highway  One  North 
Princeton,  NJ  08540 

2  Director 

National  Aeronautics  and 
Space  Administration 
Langley  Research  Center 
Langley  Station 
Hampton,  VA  25365 
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1  Director 

National  Aeronautics  and 
Space  Administration 
Manned  Spacecraft  Center 
ATTN:  Library 
Houston,  TX  77058 

1  Director 

NASA  -  Ames  Research  Center 
ATTN:  Tech  Lib 
Moffett  Field,  CA  94035 

1  Aeronautical  Research  Association 
of  Princeton,  Inc. 

50  Washington  Road 
Princeton,  NJ  08540 

1  Forrestal  Research  Center 
Aeronautical  Engineering  Lab 
Princeton  University 
ATTN:  Dr.  Hr ingen 
Princeton,  NJ  08540 

1  Northrup  Research  and  Technology  Center 
3401  W.  Broadway 
Hawthorne,  CA  90250 


1  Northrop  Research  6  Technology  Center 
ATTN :  Library 
One  Research  Park 
Palos  Verdes  Peninsula,  CA  90274 


1  .General  Electric  -  TEMPO 
ATTN:  W.  Chan 

816  State  Street 

P.0.  Drawer  QQ 

Santa  Barbara,  CA  93102 

2  Aerospace  Corporation 
ATTN:  Mr.  L.  Rubin 

Mr.  L.G.  King 

2350  E.  El  Segundo  Boulevard 
El  Segundo,  CA  90.245 
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1  Aerospace  Corporation 
ATTN:  Dr.  T.  Taylor 
P.O.  Box  92957 
Los  Angeles,  CA  90009 

1  Aircraft  Armaments  Inc. 

ATTN:  John  Hebert 

York  Road  5  Industry  Lane 
Cockeysville,  MD  21030 

2  ARES  Inc. 

ATTN:  Duane  Summers 
Phil  Conners 
Port  Clinton,  OH  43452 

1  ARO,  Inc. 

Arnold  APS, 

TN  37389 

1  BLM  Applied  Mechanics  Consultants 
ATTN:  Dr.  A.  Bores i 
3310  Willett  Drive 
Laramie,  WY  82070 

1  Boeing  Aerospace  Company 
ATTN:  Mr.  R.G.  Blaisdell 
(M.S.  40-25) 

Seattle,  WA  98124 

I  CALSPAN  Corporation 
ATTN:  E.  Fisher 
P.O.  Box  400 
Buffalo,  NY  14225 

1  Computer  Code  Consultants 
1680  Camino  Redondo 
Los  Alamos,  NM  87544 

1  Effects  Technology,  Inc. 

5383  Hollister  Avenue 
P.O.  Box  30400 

Santa  Barbara,  CA  93105 

2  Falcon  Rf(D  Company 
ATTN:  L.  Smith 

R.  Miller 

109  Inverness  Drive,  East 
Englewood,  CO  80112 


1  Falcon  R$D  Company 
Thor  Facility 
696  Fairmont  Avenue 
Baltimore, 

MD  21204 


1  FMC  Corporation  ■ 

Ordnance  Engineering  Division 
San  Jose,  CA  95114 

2  General  Electric  Company 
ATTN:  H.J.  West 

J.  Pate 

100  Plastics  Avenue 
Pittsfield,  MA  01203 

1  General  Electric  Company 
ATTN :  H . T .  Nagamat su 
1046  Cornelius  Avenue 
Schenectady,  NY  12309 

1  Kaman  -  TEMPO 
719  Shamrock  Rd. 

ATTN:  E.  Bryant 
Bel  Air,  MD  21014 

2  General  Electric  Company 

ATTN:  Armament  Systems  Department 
David  A.  Graham 
Lakeside  Avenue 
Burlington,  VT  05402 

1  President 

General  Research  Corporation 
ATTN:  Lib 
McLean,  VA  22101 

1  Goodyear  Aerospace  Corporation 
1210  Massillon  Road 
Akron,  OH  -<4315 

1  J.D.  Haltiwanger 

Consulting  Services 
B106a  Civil  Engineering  Building 
208  N.  Romine  Street 
Urbana,  IL  61801 
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1  Hercules  Inc. 

Industrial  Systems  Department 
P.0.  Box  54S 
McGregor,  TX  76657 

3  Honeywell 

Government  §  Aerospace  Froducts 
Division 

ATTN:  Mr.  J.  Blackburn 
Dr.  G.  Johnson 
Mr.  R.  Simpson 
600  Second  Street,  NE 
Hopkins,  UN  55343 

1  Hughes  Aircraft  Co, 

Bldg.  6,  MSE-125 
Centinela  6  Tenle  Streets 
Culver  City,  CA  90230 

1  Kaman-Nuclear 

ATTN:  Dr.  P.  Snow 

1500  Garden  of  the  Gods  Road 

Colorado  Springs,  CO  80907 

1  Lockheed  Missiles  §  Space  Co,  Inc. 
ATTN:  Dr.  C.E.  Vivian 
Sunnyvale,  CA  94086 

1  Lockheed  Huntsville 
P.0.  Box  U03 
Huntsville, 

AL  35809 

1  Martin  Marietta  Corp. 

Orlando  Division 
P.0.  Box  5837 
Orlando, 

FL  32805 

1  McDonnell  Douglas  Astronautics 
ATTN:  Mail  Station  21-2 
Dr.  J.  Wall 
5301  Bolsa  Avenue 
Huntington  Beach,  CA  92647 


1  Pacific  Technical  Corporation 
ATTN:  Dr.  P.K.  Feldman 
460  Ward  Drive 
Santa  Barbara,  CA  931  OS 

1  Northrup  Norair 
Aircraft  Division 
3901  W.  Broadway 
Hawthorne,  CA  ,90250 

2  Physics  International  Company 
ATTN:  Dr.  D.  Orphal 

Dr.  E.T.  Moore 
2700  Merced  Street 
San  Leandro,  CA  94577 

1  Rockwell  International 

Autometics  Missile  Systems  Division 
ATTN:  Dr.  M.  Chawls 
4300  E.  5th  Avenue 
Columbus,  OH  43216 

1  R5D  Associates 

P.0.  Box  9695 

Marina  Del  Rey,  CA  90291 

1  Science  Applications,  Inc. 

ATTN:  G.  Burghart 

201  W.  Dyer  Road  (Unit  B) 

Santa  Ana,  CA  92707 

1  Science  Applications,  Inc. 

101  Continental  Bldg.  Suite  310 
El  Segundo,  CA  90245 

1  Science  Applications,  Inc. 

2450  Washington  Avenue,  Suite  120 
San  Leandro,  CA  94577 

1  Science  Applications,  Inc. 

1710  Goodridge  Drive 
P.0.  Box  1303 
McLean,  VA  22102 


1  Northrup  Norair 
Aircraft  Division 
3901  W.  Broadway 
Hawthorne,  CA  90250 
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1  Science  Applications,  Inc. 
ATTN:  Dr.  Trivelpiece 
1250  Prospect  Plaza 

La  Jolla,  CA  92037 

2  Systems,  Science  §  Software 
ATTN:  Dr.  R.  Sedgwick 

Ms.  L.  Hageman 
P.O.  Box  1620 
La  Jolla,  CA  92037 

1  Teledyne  Drown  Engineering 
ATTN:  Mr.  John  H,  Hennings 
Cummings  Research  Park 
Huntsville,  AL  3S807 

1  SfiD  Dynamics,  Inc. 

75S  New  York  Avenue 
Huntington,  NY  11743 

1  Southwest  Research  Institute 
ATTN:  P.  Cox 
8500  Culebra  Road 
San  Antonio,  TX  78228 

1  Southwest  Research  Institute 
Lire  Research  Station 
8500  Culebra  Rd. 

San  Antonio, 

TX  78228 


2  Sout  .west  Research  Institute 
Department  of  Mechanical  Sciences 
ATTN:  Dr.  U.  Lindholm 

Dr.  IV.  Baker 
8500  Culebra  Road 
San  Antonio, 

TX  78228 

3  SRI  INtei'national 

ATTN:  Dr.  L,  Seaman 

Dr.  D.  Curran 

Dr.  D.  Shockey 

333  Ravenwood  Avenue 

Menlo  Park,  CA  94025 
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2  University  of  Arizona 

Civil  Engineering  Department 
ATTN:  Dr.  D.A.  DaDeppo 
Dr,  R.  Richard 
Tucson,  A 2  85721 

1  Brigham  Young  University 

Department  of  Chemical  Engineering 
ATTN:  Dr.  M.  Beckstead 
Provo,  UT  84601 

1  University  of  California 

Lawrence  Livermore  Laboratory 
ATTN:  Dr.  IVm,  J.  Singleton,  L-9 
P.O.  Box  808 
Livermore,  CA  94550 


2  University  of  Wisconsin-Madison 
Mathematics  Research  Center 
ATTN:  Dr.  John  Nohel 
610  Walnut  Street 
Madison, 

WI  53706 


1  University  of  Dayton 

Research  Institute 
A1TN:  S.J.  Bless 
Dayton,  OH  45406, 

2  University  of  Delaware 

ATTN:  Prof.  J.  Vinson 

Dean  J.  Greenfield 

Newark, 

DE  19711 


1  University  of  Denver 

Denver  Research  Institute 
ATTN :  Mr.  R.F.  Recht 
2390  South  University  Boulevard 
Denver,  CO  80210 

1  Drcxel  University 

Dept,  of  Mechanical  Engineering 
ATTN :  Dr .  P . C .  Chou 
32nd  and  Chestnut  Streets 
Philadelphia,  PA  19104 
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USiiR  EVALUATION  01-  REPORT 


Please  take  a  few  minutes  to  answer  the  questions  below;  tear  out 
this  sheet,  fold  as  indicated,  staple  or  tape  closed,  and  place 
in  the  mail.  Your  comments  will  provide  us  with  information  for 
improving  future  reports. 

1 .  !3RL  Report  Number _ _ 

2.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related 
project,  or  other  area  of  interest  for  which  report  will  be  used.) 


3.  How,  specifically,  is  the  report  being  used?  (Information 
source,  design  data  or  procedure,  management  procedure,  source  of 
ideas,  etc.)  _ 


4.  Has  the  .information  in  this  r.vi rl  led  to  any  quantitative 
savings  as  far  ns  mar.  -hours/ con  true  t  dollars  saved,  operating  costs 
avoided,  efficiencies  achieved,  etc.?  If  so,  please  elaborate. 


5.  General  Comments  (Indicate  what  you  think  should  be  changed  to 
make  this  report  and  future  reports  of  this  type  move  responsive 
to  your  needs,  move  usable,  improve  readability,  etc.) _ 


6,  If  you  would  like  to  be  contacted  by  the  personnel  who  prepared 
this  report  to  raise  specific  questions  or  discuss  the  topic, 
please  fill  in  the  following  information. 


Te  1  ephone  Number :  __ 
Organ! zat ion  Addrcs s : 


